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Abstract 

The maximum cardinality and maximum weight matching problems can be solved in time 0{m^/n), 
a bound that has resisted improvement despite decades of research. (Here m and n are the number of 
edges and vertices.) In this article we demonstrate that this "my^n barrier" is extremely fragile, in the 
following sense. For any e > 0, we give an algorithm that computes a (1 — e)-approximate maximum 
weight matching in 0(me~^ loge~^) time, that is, optimal linear time for any fixed e. Our algorithm is 
dramatically simpler than the best exact maximum weight matching algorithms on general graphs and 
should be appealing in all applications that can tolerate a negligible relative error. 

Our second contribution is a new exact maximum weight matching algorithm for integer-weighted bi- 
partite graphs that runs in time 0{m^/nlog N). This improves on the 0(A''my'n)-time and 0{m^/nlog{nN))- 
time algorithms known since the mid 1980s, for 1 <C logA'^ <^ logn. Here A'^ is the maximum integer 
edge weight. 

1 Introduction 

Graph matching is one of the most well studied problems in combinatorial optimization. The original mo- 
tivations of the problem were minimizing transportation costs [62l El] and optimally assigning personnel to 
job positions |29l 1107) . Over the years matching algorithms have found applications in scheduling, approxi- 
mation algorithms, network switching, and as key subroutines in other optimization algorithms, for example, 
undirected shortest paths [55] , planar max cut [221 , Chinese postman tours [35J [5T] , and metric traveling 
salesman [15j . In most practical applications it is not critical that the algorithm produce an exactly optimum 
solution. In this article we explore the extent to which this freedom — not demanding exact solutions — allows 
us to design simpler and more efficient algorithms. 

In order to discuss prior work with precision we must introduce some notation and terminology. The 
input is a weighted graph G — (V, E, w) where n = \V\ and m — \E\ are the number of vertices and edges and 
w is the edge weight function. If w assigns integer (rather than real) weights, let N be the largest magnitude 
of a weight. An unweighted graph is one for which w{e) — 1 for all e € E. A matching is a set of vertex- 
disjoint edges and a perfect matching is one in which all vertices are matched. The weight of a matching is 
the sum of its edge weights. We use mwm (and mwpm) to denote the problem of finding a maximum weight 
(perfect) matching, as well as the matching itself. We use mcm and mcpm for the cardinality (unweighted) 
versions of these problems. The mwpm problem on bipartite graphs is often called the assignment problem. 

The MWPM and mwm problems are reducible to each other. Given an instance G of mwm, let G' consist 
of two copies of G with zero-weight edges connecting copies of the same vertex. Clearly a mwpm in G' 

*This work is supported by NSF CAREER grant no. CCF-0746673 and a grant from the US-Israel Binational Science Foun- 
dation. H.-H. Su is supported by a Taiwan (R.O.C.) Ministry of Education Fellowship. Authors' emails: duanr02@gmail.com, 
pettie@umich.edu, hsinhao@umich.edu. 



1 



Table 1: Cardinality Matching 



Year 


Authors 


Time Bound & Notes 




folklore / trivial 


TflTi DLrAtxL i i hj 


1965 


Edmonds 


poly(n) 


lyoo 


vvitzgaii 06 ZjciLin 




lyoy 


Balinski 




1 Q7/1 

ly 1 D 


Kameda Munro 
Gabow 


mn or mna{m,n) or 


1976 


Lawler 




ly ( D 


Karzanov 




1 071 

ly ( 1 
1973 


Hopcroft & Karp 
Dinic & Karzanov 


VHy/n BIPARTITE 


lyou 


Micali Sz Vazirani 




1 QQ1 

lyy 1 


Gabow & Xarjan 


lyoi 


Ibarra Sz Moran 


71 L-AKUliN ALi 1 Y UiNLiY ,KAiN UUMiZr!jJJ,riiFAKi i 1 hj 


1989 


Rabin & Vazirani 


n"^ CARDINALITY ONLY,RANDOMIZED 
il HAN UUMiZjrjD 


1991 


Alt, Blum, Mehlhorn & Paul 


Uy/nm/ \ogn BIPARTITE 


1991 
1997 


Feder & Motwani 
Goldberg & Kennedy 


mJniK BIPARTITE, K = -, — ^"^o') , 


1996 


Cheriyan & Mehlhorn 


^2 _|_ 1^ BIPARTITE, w = machine word size 


2004 


Goldberg & Karzanov 




2004 
2006 


Mucha & Sankowski 
Harvey 


RANDOMIZED 



Note: Here cj < 2.376 is the exponent oi n x n matrix multiplication. 



corresponds to a pair of mwms in G. In the reverse direction, if G is an instance of mwpm with weight 
function w, find the MWM of G using the weight function w'(e) = w{e) + nN. Maximum weight matchings 
with respect to w' necessarily have maximum cardinality. Call a matching 5 -approximate, where 5 € [0, 1], 
if its weight is at least a factor 5 of the optimum matching. Let (5-MWM (and (5-MCM) be the problem of 
finding ^-approximate maximum weight (cardinality) matching, as well as the matching itself. 

Tables [T] [2] and |3] give an at-a-glance history of exact matching algorithms. Algorithms are dated 
according to their initial publication, and are included either because they establish a new time bound, or 
employ a noteworthy technique, or are of historical interest. Table |4] gives a history of approximate MCM 
and MWM algorithms. 

1.1 Algorithms for Bipartite Graphs 

The MWM problem is expressible as the following integer linear program, where x represents the incidence 
vector of the matching. 
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maximize w{e)x(e) 

eG-E 

subject to < x{e) < 1 

x{e) is an integer 



Ve e £; 



(1) 



(2) 



It is well known that in bipartite graphs the integrality requirement ([2| is redundant, that is, the basic 
feasible solutions of the LP ([T]) are nonetheless integral. See [TOl [T5] . The dual of ^ is 



minimize ui'^J-) 

subject to y{e) > w{e) 
y{u) > 

where, by definition, y{u,v) =^ y{u) + y{v) 



(3) 



(In the MWPM problem J2e=(u «') ^(^) ~ ^ holds with equality in the primal and y{u) is unconstrained 
in the dual.) Kuhn's [76LI78I publication of the Hungarian method stimulated research on this problem from 
an algorithmic perspective, but it was not without precedent. Kuhn noted that the algorithm was latent in 
the work of Hungarian mathematicians Konig and Egervaryj^ However, the history goes back even further. 
A recently rediscovered article of Jacobi from 1865 describes a variant of the Hungarian algorithm; see }90) . 
Although Kuhn's algorithm self-evidently runs in polynomial time, this mark of efficiency was noted later: 
Munkres 'M] showed that O(n^) time is sufficient. 

Kuhn's Hungarian algorithm is sometimes described as a dual (rather than primal) algorithm, due to the 
fact that it maintains feasibility of the dual (|3| and progressively improves the primal objective ([l} by finding 
augmenting paths. Gleyzal |50l (see also [8 ) gave a primal algorithm for the assignment problem in which 
the primal is feasible (the current matching is perfect) and the dual objective is progressively improved via 
weight-augmenting cycles]^ 

The search for faster assignment algorithms began in earnest in the 1960s. Dinic and Kronrod [24" gave an 
0(n'^ )-time algorithm and Edmonds and Karp ^33^ and Tomizawa |lllj observed that assignment is reducible 
to n single-source shortest path computations on a non-negatively weighted directed graphj^ Using Fibonacci 
heaps, n executions of Dijkstra's [3T] shortest path algorithm take 0{mn -t- n^logn) time. On integer 
weighted graphs this algorithm can be implemented slightly faster, in 0(mn -I- log log n) time [57t llO^ 
or 0{mn) time (randomized) [H I110| . independent of the maximum edge weight. Gabow and Tarjan j46) . 
improving an earlier algorithm of Gabow f41J , gave a scaling algorithm for the assignment problem running 
in 0{m^Jn\og{nN)) time, which is just a log(nA^) factor slower than the fastest mcm algorithm j65jj^ For 
reasonably sparse graphs Gabow and Tarjan's |46j assignment algorithm remains unimproved. However, 
faster algorithms have been developed when N is small or the graph is dense [TU [751 1103| . Of particular 
interest is Sankowski's algorithm |103j . which solves mwpm in 0{Nn'^) time, where w is the exponent of 
square matrix multiplication. 



translation of Egervary's work appears in Kuhn |77| . 

^The idea of cycle canceling is usually attributed to Robinson [lOlj . Some assignment algorithms simply do not fit the 
primal/dual mold. Von Neumann 11151 , for example, gave a reduction from the assignment problem to finding the optimum 
strategy in a zero-sum game given as an n X matrix, which can be solved in polynomial time |12|. 

^It was known that the assignment problem is reducible to n shortest path computations on arbitrarily weighted graphs. 
See Ford and Fulkerson 35 , Hoffman and Markowitz |63| . and Desler and Hakimi [191 for different reductions. 

*Gabow and Tarjan's algorithm takes a Ifungarian-type approach. The same time bound has been achieved by Orlin and 
Ahuja fSTJ using the auction approach of Bertsekas and by Goldberg and Kennedy I54| using a preflow-push approach. 
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Table 2: Weighted Matching: Bipartite Graphs 





Authors 


lime _DOunQ c6 iNotes 


1946 


Easterfield 


2"poly(n) 


iyOo 


von NGumanii 




lyOO 


ivunn 




1955 


Gleyzal 


poly(n) 


lyo / 


Munkrcs 






Balinski & Gomory 




iyoy 


Dinic & Kronrod 


n 


iy ( u 

1 071 

ly ( i 


Edmonds & Karp 
Tomizawa 


_]_ SP^ = time for one SSSP computation on 

a non-negativeiy weighted graph 


1975 




TTI Tl\Ciff ,TI ft — A- TYl 1 71 

111/11/ J-lJti J / </ U/ ^ 1 1 / L / 1 L 


1983 


Gabow 


mn^/^logA^ INTEGER WEIGHTS 

1 y 1 1 t'V/ ' t ^Vi VV iVl WIN Lj 1 ^ J. JZjvj JZjjrV VV J_J±vjTlJ. ± o 


1984 


±. J. C'U.i.i.J.Cl.LX Oo A Clii 1 dlJ- 


777 77 —1— loCTl 
111/11/ 1 / 1 lUti / 1 


1988 


V_T CX LJ vj VV Oo Cli J tXil 




1992 


Orlin & Ahuja 


mv^log(nA^) integer weights 


1997 


Goldberg & Kennedy 




1996 


Cheriyan & Mehlhorn 


n5/2 log(niV) ( '°f " ) INTEGER weights 


1999 


Kao, Lam, Sung & Ting 


Nrriy/n/K MWM ONLY, INTEGER WEIGHTS 
N{n'^ + n^/'^/w) MWM ONLY, INTEGER WEIGHTS 


2004 


Mucha & Sankowski 


Nn'^ MWM ONLY, RANDOMIZED, INTEGER WEIGHTS 


2006 


Sankowski 


iVn" RANDOMIZED, INTEGER WEIGHTS 


new 


my^logN MWM ONLY, INTEGER WEIGHTS 
rn-^n log(nA^) INTEGER WEIGHTS 


Note: 


is the maximum integer edge weight, w is the machine word size, and a — \ogn/ log(n'^/m). 



The time bounds of Johnson [TOj and Fredman and Tar j an [36| reflect faster priority queues. The 
time bound of Mucha and Sankowski ^8] follows from Kao et al.'s reduction. 



1.2 Algorithms for General Graphs 



Whereas the basic solutions to (T]3) are integral on bipartite graphs, the same is not true for general graphs. 



For example, if the graph is a unit- weighted cycle with length 2fc + 1 the mwm has weight k but ([T]) achieves 
its maximum of fc + 1/2 by setting x(e) = 1/2 for all e £ E. Let Vodd be the set of all odd-size subsets 
of V. Clearly every feasible solution to the integer linear program (Tj2) also satisfies the following odd-set 
constraints. 

a^(e) < i\B\ - l)/2 VB G Vodd (4) 

eeE{B) 

Edmonds [SUIISI] proved that if we replace the integrality constraints ^ with Q, the basic solutions to 
the resulting LP are integral]^ Edmonds' algorithm mimics the structure of the Hungarian algorithm but 
the search for augmenting paths is complicated by the presence of odd-length alternating cycles and the fact 
that matched edges must be searched in both directions. Edmonds' solution is to contract blossoms as they 
are encountered. A blossom is defined inductively as an odd-length cycle alternating between matched and 



'''In the MWPM problem x{e) = 1, for all m S V, and we have the freedom to use an alternative variety of odd-set 

constraints, namely, T,e^(u.v)eE neB,v^B ^(e) > 1, VB e Vodd- 
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Table 3: Weighted Matching: General Graphs 



Year 


Authors 


Time Bound & Notes 


lyoo 


Edmonds 


poly(n) 


1 07 A 


Gabow 




1 Q7R 


Lawler 


1 Q7R 
ly ( D 


Karzanov 


3 

n -\- TTLTi log Ti 


1 078 
ly ( o 


Cunningham & Marsh 


poly(n) 


1982 


Galil, Micali & Gabow 


mn log n 


1985 


Gabow 


TliTl ' log i V INlEGER WEIGH lb 
Nm^/n MWM ONLY, INTEGER WEIGHTS 


1989 


Gabow, Galil & Spencer 


mn log log log^ n + log n d = 2 + TO/n 


1990 


Gabow 


mn + log n 


1991 


Gabow & Tar j an 


?7?Vn.lognlog(niV) integer weights 


2006 


Sankowski 


iVn" weight only, integer weights 


2012 


Huang & Kavitha 


Nmy/n/K MWM ONLY, INTEGER WEIGHTS 
iVn" MWM ONLY, RANDOMIZED, INTEGER WEIGHTS 


Note: 

K = \0 


jy IS the maximum integer edge weight, uj is the exponent ot nx n matrix multiplication, and 
gn/ log(n^/m). 



unmatched edges, whose components are either single vertices or blossoms in their own right. Blossoms are 
discussed in detail in Section [2T1 

The fastest implementation of Edmonds' algorithm, due to Gabow (42j, runs in 0{mn + n^ logrt) time, 
which matches the running time of the best bipartite mwpm algorithm [3^. Gabow and Tarjan j47] extended 
their scaling algorithm for mwpm to general graphs, achieving a running time of 0{m^/n log n \og(nN)), 
which is the fastest known algorithm for integer- weighted graphs and nearly matches the 0{m^/n) time 
bound of the best mcm algorithms [HSJ IllSj j^ As in the bipartite case, faster algorithms for mwm and 
MWPM are known when the graph is dense or N is small. Sankowski |103j noted that the weight of the 
MWPM could be computed in 0{Nn^) time; however, it remains an open problem to adapt the cardinality 
matching algorithms of [551 EO] to weighted graphs. Huang and Kavitha generalizing [72], proved 
that MWM is reducible to N MCM computations, which, by virtue of [531 IHHl ISO], implies a new bound of 
O [N ■ min{n" , m^n log (n^ I™) I log ) ■ 

1.3 Approximating Weighted Matching 

The approximate MWM problem is remarkable in that it has been studied for decades, has practical applica- 
tions, and yet, as late as 1999, essentially nothing better than the greedy algorithm was knownj^ Moreover, 
the (1 — e)-MCM problem had been solved satisfactorily in the early 1970s. Although not stated as such, the 
0(m-\/ri)-time exact MCM algorithms [5S1 [531 [ZS [53] are actually (1 — e)-MCM algorithms running in 0(me~^) 
time. These algorithms are based on three observations (i) a maximal set of shortest augmenting paths can 
be found in linear time, (ii) augmenting along such a set increases the length of the shortest augmenting 
path, and (iii) that after k rounds of such augmentations the resulting matching is a (1 — j^)-mcm. 

Preis |97j gave a linear time ^-mwm algorithm, which improves on the greedy algorithm's 0(ra\ogn) 

^Gabow and Tarjan |47J claim a running time of 0(rn-\J n log no(m, n) log(nA'')), where the a(m, n) factor comes from an 
0{ma{rn,n)) implementation of the split-findmin data structure |40| . This can be reduced to 0(m log a(m, n)) |95| . However, 
Thorup |108| noted that split-findmin can be implemented in 0(m) time on integer-weighted graphs. 

^The greedy algorithm repeatedly includes the heaviest edge in the matching and removes all incident edges. Gabow and 
Tarjan |46l 147) observed that by retaining the 0(log(n/e)) high-order bits of the edge weights, their exact scaling algorithms 
become 0(mY^)-time (1 — e)-MWM algorithms for bipartite and general graphs. 
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Table 4: Approximate Maximum Cardinality/ Weight Matching 



Year 


Authors 


Approx. Problem 


Time Bound & Notes 




1971 
1973 


Hopcroft & Karp 
Dinic & Karzanov 


(1 — e)-MCM 




BIPARTITE 


1980 
1991 


Micali & Vazirani 
Gabow &; Tar j an 


(1 — e)-MCM 


me ^ 






folklore /trivial 


i-MWM 


TO log n 


1988 


Gabow & Tar j an 


(1 — e)-MWM 


TO-y/nlog(n/e) 


BIPARTITE 


1991 


Gabow & Tar j an 


(1 — e)-MWM 


my/ n log 71 log(7i/e) 


1999 
2003 


Preis 

Drake & Hougardy 


i-MWM 


TO 


2003 


Drake & Hougardy 


— e)-MWM 


TOe~^ 


2004 


Pettie & Sanders 


(| — e)-MWM 


TO log e^^ 


2010 
2010 


Duan & Pettie 
Hanke & Hougardy 


(| — e)-MWM 


TO log n log e^^ 


2010 


Hanke & Hougardy 


(| — e)-MWM 


TO log^ n log e^^ 


new 


(1 — e)-MWM 


me ^ log e ^ 
TOe~^ log A'' 


ARBITRARY WEIGHTS 
INTEGER WEIGHTS 



Note: N is the maximum integer edge weight and e > is arbitrary. 



running time but not its approximation guarantee. Drake and Hougardy |114j presented the first linear 
time algorithm with an approximation guarantee greater than 1/2. Specifically, they gave a (| — e)-MWM 
algorithm running in 0{me^^) time, for any e > 0. The dependence on e was later improved by Pettie 
and Sanders [SS]. These algorithms are based on a weighted version of Hopcroft and Karp's [5S] argument, 
namely that any matching whose weight-augmenting paths and cycles have at least k unmatched edges is 
necessarily a (1 — ^)-mwm. Algorithms are presented in [27l [58] [59] with different time/approximation 
tradeoffs: a (| — e)-MWM algorithm running in time 0(TOlognloge~^) and a (| ~e)-MWM algorithm running 
in 0(to log^ n log e~^) time. 

1.4 New Results 

We present the first (1 — e)-MWM algorithm that significantly improves on the 0{my/n) running times of [JHl 
147] . Our algorithm runs in 0{me~^ loge~^) time on general graphs and 0(min{TOe~^ loge~^,TOe~^ log A}) 
time on integer-weighted general graphs. This is optimal for any fixed e and near-optimal as a function of e, 
given the state-of-the-art in mcm algorithms]^ Moreover, our algorithm is as simple as one could reasonably 
hope for. Its search for augmenting paths uses depth first search |37l §8] rather than the double depth first 
search of [55]. It uses no priority queues, split-findmin structures [ID], or the blossom "shells" that arise 
from Gabow and Tarjan's [17] scaling technique. 

Our second result is a new algorithm for exact mwm on bipartite graphs running in 0{my/n\og N) time, 
which improves on [41] [46] for 1 ^ log A <C logn. According to the MWPM— j-mwm reduction, this also yields 
a new 0(TOY^log(nA)) MWPM algorithm, matching the performance of (46) . However, our algorithm can be 
used to solve MWPM directly, in log(v^A) scales rather than log(nA), which might be practically significant. 
In terms of technique, our algorithm is a synthesis of the dual (Hungarian-type) approach of Gabow and 

*Note that any (1 — e)-MWM algorithm running in 0{f{e)m) time yields an exact MCM algorithm running in 0(m ■ (/(e) -|- en)) 
time, for any e. Thus, any (1 — e)-MWM algorithm running in o{mt~^) time would improve the 0{m^/n) MCM algorithms |65l 
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Tarjan |46) and the primal approach of Bahnski and Gomory [5], among others. The ^Jn factor in our 
running time arises not from the standard blocking flow- type argument [741 165j but Dilworth's lemma |22| , 
which ensures that every partial order on n elements contains a chain or anti-chain with size ^Jn. Dilworth's 
lemma has also been used in Goldberg's [52] single-source shortest path algorithm. 

1.5 Remarks on Approximate Weighted Matching and Its Applications 

Our focus is on algorithms that accept arbitrary input graphs and that give provably good worst-case 
approximations. These twin objectives are self-evidently attractive, yet nearly all work (prior to Preis ^7\) on 
approximate weighted matching focused on specialized cases or weaker approximation guarantees. Early work 
on the problem usually considered complete bipartite graphs, and confirmed the efficiency of heuristics either 
experimentally or analytically with respect to inputs over some natural distribution [HI I107( l87l \79\ \80\ [5] . 
See Avis [6 for a more detailed discussion of heuristics. 

Most work in the area considers graphs defined by metrics, often Euclidean metrics. Reingold and 
Tarjan jlOOj proved that the greedy algorithm for metric MWPMjhas an approximation ratio of « 

Goemans and Williamson [51] gave a 2-approximation for metric MWPM that can be implemented 
in O(n^) time |44j, or O(mlog^n) time [16] in metrics defined by m-edge graphs. The Euclidean mwpm 
comes in two flavors: the monochromatic version is given 2n points and the bichromatic version is given 2n 
points, n of which are colored blue, the rest red, where the matching cannot include monochromatic edges 
Varadarajan and Agarwal |112| gave (1 -I- e)-MWPM algorithms for the mono- and bichromatic variants 
running in time 0(rtpoly(e~^ logn)) and 0(n'^/^poly(e~^ logn)), respectively. Other time-approximation 
tradeoffs for the bichromatic variant are possible [Il[13], including an O(npoly(logrt)) time algorithm for 
0(l)-approximating the weight of the mwpm [55]. Some work considers the even more specialized case of 
Euclidean matching in the unit square, which allows for algorithms that guarantee absolute upper bounds 
on the weight of the matching; see [531 [Ml IS] and the references therein. 

There are several applications of MWM (on general or bipartite graphs) in which one would gladly sacrifice 
matching quality for speed. In input-queued switches packets are routed across a switch fabric from input 
to output ports. In each cycle one partial permutation can be realized. Existing algorithms for choosing 
these matchings, such as iSLIP [84j and PIM [3], guarantee ^-MCMs and it has been shown [551133] that 
(approximate) MWMs have good throughput guarantees, where edge weights are based on queue-length. See 
also [531 11061 1105j . Approximate mwm algorithms are a component in several multilevel graph clustering 
libraries!^ (PARTY, for example, builds a hierarchical clustering by iteratively finding and contracting 
approximate MWMs; see [98].) Approximate MWM algorithms are used as a heuristic preprocessing step in 
several sparse linear system solvers [CT [1041 156j . The goal is to permute the rows/columns to maximize 
the weight on or near the main diagonal. 

1.6 Organization 

Section [2] introduces some notation, states well known properties of augmenting paths and blossoms, and 
reviews Edmonds' optimality conditions for weighted matching. In Section [s] we present our (1 — e)-MWM 
algorithm and in Section [2] we present exact algorithms for bipartite mwm and mwpm. 

2 Preliminaries 

We use E{H) and V{H) to refer to the edge and vertex sets of H or the graph induced by H , that is, V{E') 
is the set of endpoints oi E' C E and E{V') is the edge set of the graph induced by V' C V. A matching M 
is a set of vertex-disjoint edges. Vertices not incident to an M edge are free. An alternating path (or cycle) 

®For metric inputs let MWPM be the minimum weight perfect matching problem. 

^"^The weight of the bichromatic MWPM is also known as the earth mover distance between the red and blue points. 
"E.g., METIS .73 , PARTY |98], PT-SCOTCH [94| CHACO [61], JOSTLE [116], and KaFFPa/KaFFPaE fg4llT02l . 
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is one whose edges alternate between M and E\M . An alternating path P is augmenting if P begins and 

ends at free vertices, that is, M ® P =^ (^\^') U (P\M) is a matching with cardinality \M ® P\ = \M\ + 1. 

When we seek (1 — e) approximate solutions, we can afford to scale and round edge weights to small 
integers. To see this, observe that the weight of the mwm is at least Wmax = max{w(e) | e € E{G)}. It 
suffices to find a (I — e/2)-MWM M with respect to the weight function i()(e) — \ w{e)/^\ where 7 = t-Wmax/n. 
Note that w{e) — 7 < 7 • 'w{e) < ■w{e) for any e. It follows from the definitions that: 



w{M) > 7 • w{M) 

>7-(I-e/2)iZ;(Ar) 

> (I-e/2)(w(M*)-77i/2) 

= {l-e/2){w{M*)~e-Wn,,^/2) 

> (I - e)w{M*) 



Defn. of w 
Defn. of M, M* is the mwm 
Defn. of w, |A/*| < n/2 
Defn. of 7 
Since w{M*) > 



Since it is better to use an exact mwm algorithm when e < 1/n [151 SZIj we assume, henceforth, that 
w : E {1,2, . . . , N}, where < is the maximum integer edge weight. For notational convenience we 
also assume that iV is a power of 2. 



2.1 Blossoms and the LP Formulation of MWM 

The dual LP of ([T]4| is 



minimize 



u<£V{G) seVodd 
subject to yz{e) > w{e) 

y{u) > 0,z{B) > 



Ve e E{G) 
Wu e V{G),yB eVodd 



where, by definition, yz{u, v) = y{u) + y{v) + z{B) 



seVodd, 

{u,v)eE{B) 

Despite the exponential number of primal constraints and dual z-variables, Edmonds demonstrated that 
an optimum matching could be found in polynomial time without maintaining information (z-values) on 
more than n/2 elements of Vodd at any given time. At intermediate stages of Edmonds' algorithm there is a 
matching M and a laminar (nested) subset D, C Vodd, where each element of Q is identified with a blossom. 
Blossoms are formed inductively as follows. If 1; g ^ then the set {v} is a trivial blossom. An odd length 
sequence {Aq,Ai, . . . , A^) forms a nontrivial blossom B = IJi^i {^i} are blossoms and there is a 

sequence of edges eo, . . . ,ei where E Ai x yl^+i (modulo ^ + 1) and e M if and only if i is odd, that is, 
Aq is incident to unmatched edges eo,e^. See Figure [l] The base of blossom B is the base of Aq; the base 
of a trivial blossom is its only vertex. The set of blossom edges Eb are {eo, . . . , } and those used in the 
formation of Aq, . . . , Ai. The set E{B) =^ E {B x B) may, of course, include many non-blossom edges. A 
short proof by induction shows that |-B| is odd and that the base of B is the only unmatched vertex in the 
subgraph induced by B. 

Matching algorithms represent a nested set ft of active blossoms by rooted trees, where leaves represent 
vertices and internal nodes represent nontrivial blossoms. A root blossom is one not contained in any other 
blossom. The children of an internal node representing B are ordered according to the odd cycle that formed 
B, where one child is distinguished as containing the base of B. As we will see, it is often possible to treat 
blossoms as if they were single vertices. Let the contracted graph G/fl be obtained by contracting all root 
blossoms and removing spurious edges. To dissolve a root blossom B means to delete its node in the blossom 
forest and, in the contracted graph, to replace B with individual vertices Aq, . . . ,Ai. Lemma [l] summarizes 
some useful properties of the contracted graph. 
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Lemma 1. Let n be a set of blossoms with respect to a matching M . 

1. If M is a matching in G then M/VL is a matching in G/VL. 

2. Every augmenting path P' relative to AI/D, in G/fl extends to an augmenting path P relative to M in 
G. (That is, P is obtained from P' by substituting for each non-trivial blossom vertex B in P' a path 
through Eb- See Figure^^a,b).) 

3. If P is an augmenting path and P/Q is also an augmenting path relative to M/fl, then O remains a valid 
set of blossoms (possibly with different bases) for the augmented matching M Q) P. See Figure^a,b). 

4-. The base u of a blossom B ^ U, uniquely determines a maximum cardinality matching of Eb, having 
size {\B\ - l)/2. See Figure\^a,b). 

Implementations of Edmonds' algorithm grow a matching M while maintaining Property[l] which controls 
the relationship between Af , f2 and the dual variables. 

Property 1. (Complementary Slackness Conditions) 

1. Nonnegativity: z{B) > for all B E Vodd and y{u) > for all u E V{G). 

2. Active Blossoms: il contains all B with z{B) > and all root blossoms B have z{B) > 0. (Non-root 
blossoms may have zero values.) 

3. Domination: yz{e) > w{e) for all e E E. 

4. Tightness: yz{e) = w{e) when e G A/ or e G Eb for some _B G 51. 

If the y- values of free vertices become zero, it follows from domination and tightness that Af is a maximum 
weight matching, as the following short proof attests. Here M* is any maximum weight matching. 



w{M) = ^ w{e) 

eSAf 

= ^ yz{e-) tightness 



E yH + E^^-^(^) Note E y(^)= E y( 

uGV(G) Ben uev{G) iiev{M) 



u 



> E + E l-^(^) ^ ^^* \ ' ^ non-negative 

u&v{M') Ben 

— E] ^ w{M*) domination 

eeM* 

3 A Scaling Algorithm for Approximate MWM 

Our algorithm maintains a dynamic relaxation of complementary slackness. In the beginning domination is 
weak but becomes progressively tighter at each scale whereas tightness is weakened at each scale, though 
not uniformly. The degree to which a matched edge or blossom edge may violate tightness depends on when 
it last entered the blossom or matching. 

Recall that N is the maximum integer edge weight. The parameter e' = 0(e) will be selected later to 
guarantee that the final matching is a (1 — e)-MWM. Henceforth, assume that > 1 and e' < 1/4 are powers 
of two. Define Sq ~ e'N and 6i = 6o/2^. At scale i we use the truncated weight function Wi{e) — 6i lw{e)/6i\ . 
Note that ^^+1(6) = Wi{e) or Wi{e) + 
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(a) (b) 

Figure 1: Thick edges are matched, thin unmatched, (a) A blossom Bi = {ui,U2, B2,us,uq,uiq, B^) with 
base ui containing non-trivial sub-blossoms B2 = ("3, U4, U5, ug, M7) with base U3 and ^3 = (un, "12, W13) 
with base uu. Vertices mi5,mi6, and un are free. The path (ttig, M2, "3, "7, U6j "Sj W4, W17) is an exam- 
ple of an augmenting path that exists in G but not G/Bi, the graph obtained by contracting Bi. (b) 
The situation after augmenting along {ui5,ui4, Bi,un) in G/Bi, which corresponds to augmenting along 
(mi5, ui4, Ml, W2, "3, "7, "61 "Sj "4j W17) in G. After augmentation Bi and B2 have their base at U4. 



Property 2. (Relaxed Complementary Slackness) There are L + 1 scales numbered 0, . . . , L, where L = 
logiV. Let i E [0,L] be the current scale. 

1. Granularity: z(B) is a nonnegative multiple of 6i, for all B e Vodd, and y{u) is a nonnegative multiple 
of (5i/2, for aU u e V{G). 

2. Active Blossoms: il contains all B with z{B) > and all root blossoms B have z{B) > 0. (Non-root 
blossoms may have zero z- values.) 

3. Near Domination: yz(e) > Wi{e) — 6i for all e € E. 

4. Near Tightness: Call a matched or blossom edge type j if it was last made a matched or blossom edge 
in scale j < i. (That is, it entered the set M U Usgo in scale j and has remained in that set, even 
as M and Q evolve as augmenting paths are found and blossoms are formed and dissolved.) If e is such 
a type j edge then yz{e) < Wi{e) + 2{Sj — Si). 

5. Free Vertex Duals: The y- values of free vertices are equal and strictly less than the ?/- values of matched 
vertices. 

Lemma [2] allows us to measure the quality of a matching M, given duals y and z satisfying Property [2] 

Lemma 2. Let M be a matching satisfying Property^ at scale i and let M* be a maximum weight matching. 
Let f be the number of free vertices, each having y-value (j), and let e be such that yz{e) — w{e) < e ■ w{e) for 
alleG M. Thenw{M) > {l + e)-^{w{M*)-25^\M*\- fcf)). Lfi = L and (j) = then M is a (l-e'-e)-MWM. 
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Proof. The claim follows from Property [2j 



(Af) = w{e) defn. of w{M) 

eeM 

> (1 + e)^^ yz{e) near tightness, defn. of e 



= (l + rM E + E • ^(^) defn. of 

>(i + rM E yN + El^(^)n^-^*l-^(^)-/'^) (5) 

\uev{M') Ben J 

> (1 + f")"' f E y^(^) - ) def^- of 

> (1 + e)"^ (w,(Af*) - f(t)~6,- \M*\) near domination 

> (1 + e)"^ (w(M*) - f4>- 25, ■ \M*\) defn. of w. 

Inequality ([5| follows from several facts, namely, that no matching can contain more than {\B\ — l)/2 edges 
in that V{M*)\V{M) contains only free vertices (with respect to A/), whose y-values are 0, and that 
y- and z-values are nonnegative. Note that the last inequality is loose by <5i|Af*| \i i = L since in that case 
Wl = w. 

The integrality of edge weights imphes that w{M*) > |Af*|. If « = i and = then 6l = = e' 

and w{M) > (1 + e)-i(w(Af*) - 5l\M*\) > (1 + €)-^{l ~ e')w{M*) > (1 - e' - i)w{M*), that is, Af is a 
(1 — e' — e)-MWM. □ 

While not suggesting an algorithm per se, Lemma[2]tells us which invariants our algorithm must maintain 
and when it may halt with a (1 — e)-MWM. For example, if at the last scale L we have 6l < e/2 and for any 
type j edge e € Af, 6j < (e/4)w(e) (i.e., e — e/2) then as soon as y-values at free vertices reach zero, the 
current matching must be a (1 — e)-MWM. 



3.1 The Scaling Algorithm 

Initially Af = 0, = 0, and y{u) = N/2 — 6q/2 for all u £ V, which clearly satisfies Property [2] for scale i — 0, 
since yz(e) ~ 2{N/2—5q/2) > wo{e) — 5Q. The algorithm, given in Figure [2] consists of scales 0, L = log iV, 
where the purpose of scale i is to halve the y-values of free vertices while maintaining Property [2j In each 
iteration of scale i the algorithm (1) augments a maximal set of augmenting paths of eligible edges, (2) finds 
and contracts blossoms of eligible edges, (3) performs dual adjustments on y- and z-values, and (4) dissolves 
previously contracted root blossoms if their z-values become zero. Each dual adjustment step decrements 
by Si/2 the y-values of free vertices. Thus, there are roughly {N/2^)/{Si/2) = 0{e^^) iterations per scale, 
independent of i. The efficiency and correctness of the algorithm depend on eligibility being defined properly. 

Definition 1. At scale i, an edge e is eligible if at least one of the following hold: 

(i) e e Eb for some B £ fl. 

(ii) e ^ Af and yz(e) = Wi(e) — Si. 

(iii) e G M and yz(e) — Wi{e) is a nonnegative integer multiple of Si. 

Let Eeiig be the set of eligible edges and let Gehg — {V, Eeiig)/fl be the unweighted graph obtained by 
discarding ineligible edges and contracting root blossoms. 
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Criterion (i) for eligibility simply ensures that an augmenting path in Geiig extends to an augmenting 
path of eligible edges in G. A key implication of Criteria (ii) and (iii) is that if P is an augmenting path in 
Gehg, every edge in P becomes ineligible in (AZ/fi) ® P. This follows from the fact that unmatched edges 
must have yz{e) — Wi{e) < whereas matched edges must have yz{e) — Wi{e) > 0. Regarding Criterion (iii), 
note that Property [2] (granularity and near domination) implies that yz{e) — Wi{e) is at least —5i and an 
integer multiple of 5i/2. 



3.2 Analysis and Correctness 

Lemma [5] states that the algorithm maintains Property [2] after each of the 0{e~^) Dual Adjustment steps in 
each scale. Lemma [3] establishes the critical fact that all augmenting paths and paths from free vertices to 
blossoms are eliminated from Geiig before each Dual Adjustment step. 

Lemma 3. After the Augmentation and Blossom Shrinking steps Geiig contains no augmenting path, nor 
is there a path from a free vertex to a blossom. 

Proof. Suppose there is an augmenting path P in G^^ug after augmenting along paths in Since is 
maximal, P must intersect some P' e 5* at a vertex v. However, after the Augmentation step every edge in 
P' will become ineligible, so the matching edge {v,v') S M is no longer in Geiig, contradicting the fact that 
P consists of eligible edges. Since fl' is maximal there can be no blossom reachable from a free vertex in 
Gehg after the Blossom Shrinking step. □ 

Lemma |4] guarantees that all y-values updated in a Dual Adjustment step have the same parity as a 
multiple of Si/2. In the proof of Lemma [S] this fact is used to argue that if both endpoints of an edge e have 
their j/-values decremented, then yz{e) is a multiple of Si. 

Lemma 4. Let R C V{Geiig) be the set of vertices reachable from free vertices by eligible alternating paths, 
at any point in scale i. Let R C V{G) be the set of original vertices represented by those in R. Then the 
y-values of R-vertices have the same parity, as a multiple of Si/2. 

Proof. Assume, inductively, that before the Blossom Shrinking step, all vertices in a common blossom have 
the same parity, as a multiple of Si/2. Consider an eligible path P = {Bq, Bi, . . . , B^.) in G^ug, where the 
{Bj} are either vertices or blossoms in and Bo is unmatched in Geiig. Let (wq, f 1), (ui, W2), ■ ■ • , {uk-i, Vk) 
be the G-edges corresponding to P, where Uj,Vj € Bj. By the inductive hypothesis, Uj and vj have the 
same parity, and whether {uj,Vj^i) is matched or unmatched, Definition [l] stipulates that yz{uj,Vj+i)/Si is 
an integer, which implies y{uj) and y(uj+i) have the same parity as a multiple of Si/ 2. Thus, the y-values of 
all vertices in i?o U • • • U B^. have the same parity as a free vertex in Bq, whose y- value is equal to every other 
free vertex, by Property [2j[5]) . Since new blossoms are formed by eligible edges, the inductive hypothesis is 
preserved after the Blossom Shrinking step. It is also preserved after the Dual Adjustment step since the 
y-values of vertices in a common blossom are incremented or decremented in lockstep. This concludes the 
induction. □ 

Lemma 5. The algorithm preserves Property^ 



Proof. Property [2|[5]) (free vertex duals) is obviously maintained as only free vertices have their y-values 
decremented in each Dual Adjustment step. Property [2|[2]) (active blossoms) is also maintained since all the 
new root blossoms discovered in the Blossom Shrinking step are contained in Vout and will have positive z- 
values after adjustment. Furthermore, each root blossom whose z-value drops to zero is dissolved, after Dual 
Adjustment. At the beginning of scale i all y- and z-values are integer multiples of Si/2 and Si, respectively, 
satisfying Property [2|T]) (granularity). This property is clearly maintained in each Dual Adjustment step. If 
e ^ M is placed in M during an Augmentation step or placed in Usen during a Blossom Shrinking step 
then e is type i and yz{e) = Wi{e) — Si < Wi{e), which satisfies Property [2j[4]) . 

It remains to show that the algorithm maintains Property [2 13|4 1 (near domination and near tightness). 



First consider the dual adjustments made at the end of scale i (the last line of pseudocode in Figure [2]) 
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Initialization: 



M ^ 
f7 ^ 
5o ^ e'N 



no matched edges 
no blossoms 

e' = 6(e); w.l.o.g., e' ,N are powers of 2 
satisfies Property [2j[3]) 




Execute scales i — . . . , L = log N and return the matching M. 



Scale »: 



— Repeat the following steps until y-values of free vertices reach iV/2*+^ — Si/2, if i € [0, L), 
or until they reach zero, if i = L. 

(1) Augmentation: 

Find a maximal set 5* of augmenting paths in GeHg and set M ^ M G) (Upe*^)- 
Update Geiig- 

(2) Blossom Shrinking: 

Let Vout C V{Geiig) be the vertices (that is, root blossoms) reachable from free vertices 
by even-length alternating paths; let fi' be a maximal set of (nested) blossoms on Vout- 
(That is, if {u,v) € E{GeUg)\M and u,v € Vouti then u and v must be in a common 
blossom.) Let Vin C V{Geiig)\Vout be those non-Vout-vertices reachable from free 
vertices by odd-length alternating paths. Set z{B) 4- for i? G f2' and set 57 <— U D! . 
Update Geiig- 

(3) Dual Adjustment: 

Let Vin , Vout C be original vertices represented by vertices in Vin and Vout ■ The y- 
and z-values for some vertices and root blossoms are adjusted: 



(4) After dual adjustments some root blossoms may have zero z-values. Dissolve such 
blossoms (remove them from $7) as long as they exist. Note that non-root blossoms are 
allowed to have zero z-values. Update Gf-ug by the new fi. 

— Prepare for the next scale, if i € [0, L): 



y{u) ^ y{u) - 5i/2, for all u e Vout- 
y{u) -s- y{u) + 5i/2, for all u e Vn- 

z{B) -s— z{B) + Si, if S G 17 is a root blossom with B C Vout- 
z{B) z{B) — Si, if S e 57 is a root blossom with B C Vin- 



Si+i ^ Si/2 

y{u) -s- y{u) + Si+i, for all u e V{G). 



Figure 2: A (1 — e)-approximate mwm algorithm 
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Let e — {u,v) be an arbitrary edge and let yz and yz' be the function before and after dual adjustment. It 
follows that 

yz'{e) = yz{e) + 25.,+i y{u), y{v) incremented by S.^+i 

> Wi{e) — Si + 2Si-^.l near domination at scale i 

> Wi+i{e) ~ Si+i Wi{e) > Wi+i{e) - Si+i 

That is, Property [2|[3]) is preserved. If e e AI (J Usesi is a type j edge, then at the end of scale i 
Property [2||4]) is also preserved since 

yz'{e) = yz{e) + 2(5,+i < w,(e) + 2{6j - S^) + 2S^+i < w^+i{e) + 2{dj - S,+i) 

The first inequality follows from Property [2||4| at scale i and the second inequality from the fact that 
Wi(e) < Wi+i(e) and (5i = 2(5i+i. 

Now consider a Dual Adjustment step. If neither u nor v is in Vin U Vout or if u, v are in the same root 
blossom in fl, then yz{e) is unchanged, preserving Property[2] The remaining cases depend on whether (u, v) 
is in M or not, whether (u, v) is eligible or not, and whether both u,v ^ Vin U Vout or not. 

Case 1: e ^ M, u,v £ Vn U Vout If e is ineligible then yz{e) > Wi{e) — 5i. However, by Lemma [i] (parity of 
2/-values) we know {yz{e) — Wi{e)) / Si is an integer, so yz{e) > Wi{e) before adjustment and yz{e) > Wi{e) — Si 
afterward (which could occur if both u,v G Vout), thereby preserving Property [ijis]) . If e is eligible then at 
least one of u,v is in Vim otherwise another blossom or augmenting path would have been formed, so yz{e) 
cannot be reduced, which also preserves Property [2|[3]) . 

Case 2: e e M, u,v E Vn U Vout Since u,v G Vn U Vout, Lemma [i] (parity of y- values) guarantees that 
{yz(e) — 'Wi{e))/Si is an integer. The only way e can be ineligible is if ?/2;(e) — Wi{e) — Si andu,?; G Vin, hence 
yz{e) — Wi{e) after dual adjustment, which preserves Property [2| |3|4[ ). On the other hand, if e is eligible then 
u e Vin and v G Vout- It cannot be that u,v G Vout, otherwise e would have been included in an augmenting 
path or root blossom. In this case yz{e) is unchanged, preserving Property [2| |3|4[ ). 

Case 3: e ^ M, v ^ Vn U Vout If e is eligible then u e Vn and yz{e) will increase. If it is ineligible then 
yz(e) > Wi{e) — Si/2 before adjustment and yz(e) > Wi{e) — Si afterward. In both cases Property [2j|3]) is 
preserved. 



Case 4: e €E M, v ^ Vin U Vout It must be that e is ineligible, so u G Vm and yz{e) — Wi{e) is either 
negative or an odd multiple of 5^/2. If e is type j then, by Property [2{ Tp ) (granularity and near tightness), 
yz[e) < Wi{e) + 2{5j — Si) — Si/2 before adjustment and yz{e) < Wi{e) + 2{Sj — Si) afterward, preserving 
Property!^. □ 



Recall that Lemma [2] stated that the final matching will be a (1 — 0(e))-MWM if Sl = 0{e), free vertices 
have zero j/- values, and yz{e) — w{e) = 0{e) ■ w(e). Lemmas [6] and [t] establish these bounds. 

Lemma 6. Let i < L be the scale index. Then 

1. For i < L, all edges eligible at any time in scales through i have weight at least N/2'''^^ + Si. 

2. For any i, if e € M then yz{e) < (1 + 4e')w(e). 
Proof. We prove the parts separately. 



Part 1 The last search for augmenting paths in scale i begins when the y- values of free vertices are iV/2'+^, 
and strictly less than y- values of other vertices, by Property [2j[5| . An unmatched edge e = {u,v) can only 
be eligible at this scale if yz{e) — Wi{e) ~ Si < w{e) — Si. Hence w{e) > y{u) + y{v) + Si > N/2^'^^ + Si. 
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Part 2 Let e be a type j edge in M during scale i. Property [2||4|) states that yz{e) — Wi{e) < 2{Sj — 5i). 
Since Wi{e) < w{e) it also follows that yz{e) — w{e) < 26j — 2Si < 2Sj — e'N/2^~^. By part 1, a type j edge 
must have weight at least N/2^^^ + Sj, hence yz{e) — ■w{e) < 4e' • w{e). □ 

Lemma 7. After scale L — logN, M is a {1 ~ 5e')-MWM. 

Proof. The final scale ends when free vertices have zero y- values. Property [2j[3]) holds with respect to 6l = 
(5o/2^ = e'N/2^ = e' and Lemma [g] states that yz{e) < (1 + 4e')w(e). By Lemma [2) w{M) > (1 - 
5e')w{M*). □ 

Theorem 1. j4 (1 — e)-MWM can be computed in time 0{me~^ logiV). 

Proof. Each Augmentation and Blossom Shrinking step takes 0{m) time using a modified depth-first search [Tfl 
§8]. (Finding a maximal set of augmenting paths is significantly simpler, conceptually, than finding a maxi- 
mal set of minimum- length augmenting paths, as is done in ^86 } I113j .) Each Dual Adjustment step clearly 
takes linear time, as does the dual adjustment at the end of each scale. Scale i < L = logA^ begins with 
free vertices' y- values at iV/2'+^ — 5i and ends with them at iV/2'+^ — Si. Since y- values are decremented 
by 6i/2 in each Dual Adjustment step there are exactly (iV/2'+^)/((5i/2) = N/{2So) = €'~^/2 such steps. 
The final scale begins with free vertices' y-values at iV/2^+^ — Sl and ends with them at zero, so there are 
fewer than {N/2^'^^)/{Sl/2) = e'~^ Dual Adjustment steps. Lemma guarantees that the final matching is 
a (1 — e)-MWM for e' < e/5. Hence, the total running time is 0(me~MogAf). □ 

3.3 A Linear Time Algorithm 

Our 0(TOe^^ log A^)-time algorithm requires few modifications to run in linear time, independent of N. In 
fact, the algorithm as it appears in Figure [2] requires no modifications at all: we only need to change the 
definition of eligibility and, in each scale, refrain from scanning edges that cannot possibly be eligible or part 
of augmenting paths or blossoms. In light of Lemma |6j[I]) it is helpful to index edges according to the first 
scale in which they may be eligible. 

Definition 2. Define /ij — iV/2'+^ + 5,;, for i < L, and /i^ — 0. Define scale(e) to be the i such that 
w{e) e 

Definition [3] redefines eligibility. The differences with Definition [T] are underlined. 
Definition 3. At scale i, an edge e is eligible if at least one of the following hold: 

1. e e Eb for some B <E Q. 

2. e ^ M and yz{e) ~ Wi(e) — Si. 

3. e G M and yz{e) — Wi(e) is a nonnegative intelligent multiple of Si. Furthermore, scale(e) > « — 7, 
where 7 '= loge'^^. 

Let Eeiig be the set of eligible edges and let GeUg = {V,EeUg)/^ be the unweighted graph obtained by 
discarding ineligible edges and contracting root blossoms. 



Lemma 8. Using Definition^ of eligibility rather than Definition^ Property 2[ is maintained and 



Property^^ (near tightness) holds in the following weaker form. Let e € MUlJ^en be a type j edge with 
scale(e) — i. Then yz{e) < Wk{e)+2{Sj —Sk) at any scale k G [i, i + ^] and yz(e) < Wk{e) + SSi < (l-|-6e')w(e) 
for k > i + "f. 



Proof. In scales i through i + 7 Property [2||4|) is maintained as the two definitions of eligibility are the 
same. At the beginning of scale t = i + 7 + l, eisno longer eligible and the y- values of free vertices are 
N /2*-^^ — St/2. From this moment on, the y-values of free vertices are incremented by a total of X]/>t+i 
(the dual adjustments following scales t through L — 1) and decremented a total of A^/2*+^ ^'^t/2 + X]/>t+i 
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(in the Dual Adjustment steps following searches for augmenting paths and blossoms). Each adjustment to 
y-values by some quantity A may cause yz{e) to increase by 2A. This clearly occurs in the dual adjustments 
following each scale as y{u) and y{v) are incremented by A. Following a search for blossoms it may be that 
u,v G Vin, which would also cause y{u) and y{v) to each be incremented by A. Note that y{u), y{v) cannot 
be decremented in scales t through L; if either were in Vout after a search for blossoms then e would have 
been eligible, a contradiction. Hence Property [2j|3]) (near domination) is maintained for e. Putting this all 
together, it follows that at any scale k>t — i + j + l = i + log e'^^ + 1, 



yz{e) < Wk{e) + 2(6^ - 4) + 2 • N/2'+^ - St/2 




< Wk{e) + 25, + 2 {5i+2 + |<5t) j > i, defn. of t 
= Wk{e) + (2 + 1/2 + ?.e'/2)5, defn. of 5t = e'5,/2 
<Wk{e) + i5, e'<l/3 

< (1 + 6e')w(e) w{e) > Wkie) > N/2'+^ = (2e')"^ • 5i 

which proves the claim. Note that the inequality yz{e) < Wk{e) + 'iSi also holds in scales A: g [4,1 + 7]. If ^ 
is type j >i then yz{e) < Wk{e) + 2{5j — 5k) < Wfc(e) + 25i. This fact will be used in Lemma[9[ □ 

The algorithm will deliberately ignore unmatched edges that may still be eligible according to Definition|3] 
This will be justified on the grounds that such edges must be adjacent to matched ineligible edges, and 
therefore cannot be contained in an eligible augmenting path. Lemma |9] will be used to argue that scale-i 
edges can be safely ignored after scale z + 7 + 2. 

Lemma 9. Let ei = {u,v) be an edge with scale(ei) = i and let eo = {u',u) and ei — {v,v') be the M- 
edges incident to u and v at some time after scale i. Then at least one of eo and 62 exists, say cq, and 
scale(eo) < i + 2. 

Proof. Following the last Dual Adjustment step in scale i the y- values of free vertices are A^/2*+^ — 5i/2. It 
cannot be that both u and v are free at this time, otherwise yz{ei) = y{u) + y{v) = iV/2'+^ — 5i = fXi — 25i < 
Wi{ei) — Si, violating Property [2||3]) (near domination). Hence, either u or t; is matched for the remainder of 
the computation. If ei is matched the claim is trivial, so, assuming the claim is false, whenever 60,62 exist 
we have scale(eo), scale(62) > i + 3. That is, w(6o), ^(62) < fJ.i+2 = N/2^^^ + (5^+2- 

It cannot be that 61 is in a blossom without eo or 62 also being in the blossom. For any I e {0, 1, 2} let 
C J7 be the blossoms containing ei at a given time. The laminarity of blossoms ensures that either Bi C Bq 
or Bi C B2- Suppose it is the former, that is, cq exists and 62 may or may not exist. Then, if the current 
scale is fc > i + 3, by Property [2|j3]) (near domination) yz{ei) — y{u) + y{v) + X^seSi ^{^) ^ Wfe(ei) — 5k- By 
Lemma[8]y(u) + J^BeBi ^(^) < y^(eo) < Wfc(eo) + 3(5^+3 and, if 62 exists, y{v) < yz{e2) < Wfc(e2) + 3(5^+3. 
These inequalities follow from the definition of yz, the containment Bi C Bq and the fact that eo and 62 can 
only be at scale i + 3 or higher. Without loss of generality we can assume y{u) + J2bgBi ^i^) — ui'^)- (If ^2 
exists and y{v) > y{u) + X^seei -^(^) ^^^^^ ^2 takes the role of eo below.) Putting these inequalities together 
we have 



Wk 



(ei) < y{u) + y{v) + ^ z{B) + 4 



BeBi 



<2iy{u)+ ^(^)) 
\ BeBi ) 

<2(y0(eo)) + 4 

< 2(M;fc(eo) + 3(5,+3) + 4 

< 2w(eo) + 7^4+3 



near domination 



SlCSo 
Lemma |8] 
/c > i + 3, e' < 1/3 
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and therefore 



w{eo) > ^{wkiei) - Si) TJ^+s < Si 

> N/2'+^ scale(ei) = i, Wfc(ei) > m ^ N/2'+^ + S^ 

> 7V/2^+3 + S,+2 = ^l^+2 

This contradicts the fact that scale(eo) > i + 3, since, by definition, such edges have u'(eo) < ^ii+2- D 
Theorem 2. A {1 — e)-MWM can he computed in time 0{me^^ loge^^). 

Proof. We execute the algorithm from Figure [2] where Geiig refers to the ehgible subgraph as defined in 
Definition [sj We need to prove several claims: (i) the algorithm does, in fact return a (1 — e)-MWM for 
suitably chosen e' = 0(e), (ii) the number of scales in which an edge could conceivably participate in an 
augmenting path or blossom is loge"""^ + 0(1), and (iii) it is possible in linear time to compute the scales in 
which each edge must participate. Part (i) follows from Lemmas [2] and [sj Since yz{e) < (1 + 6e')w{e) for 
any e € M and Sl — e', Lemma [2] implies that M is a (1 — e)-MWM when e' < e/7. 

Turning to part (ii), consider an edge e with scale (e) — i. By Lemma [6|T]) e can be ignored in scales 
through i—1. If e = {u, v) £ M then, according to Definition|3j e will be ineligible in scales i + 7 + 1 through 
log N. After scale i + 7 no augmenting path or blossom can contain e, so we can commit it to the final 
matching and remove from consideration all edges incident to u or v. Now suppose that e ^ M at the end 
of scale i + 7 + 2. Lemma [o] states that either m or w is incident to a matched edge eo with scale(eo) < « + 2, 
which by the argument above, will be committed to the final matching, thereby removing e from further 
consideration. Thus, to faithfully execute the algorithm we only need to consider e in scales scale(e) through 
scale(e) +7 + 2, that is, 7 + 3 = log e'^^ + 3 < log e^^ + 6 scales in total. 

We have narrowed our problem to that of computing scale(e) for all e. This is tantamount to computing 
the most significant bit (MSB(a;) = [log2 x\) in the binary representation of w{e). Once the MSB is known, 
scale(e) can be just one of two possible values. MSBs can be computed in a number of ways using standard 
instructions. It is trivial to extract MSB(a;) after converting x to floating point representation. Fredman 
and Willard ^\ gave an 0(1) time algorithm using unit time multiplication. However, we do not need to 
rely on floating point conversion or multiplication. In Section [2] we showed that without loss of generality 
log TV < 21ogn. Using a negligible 0{n^) space and preprocessing time we can tabulate the answers on 
(3 ■ logn-bit integers, where /? < 1, then compute MSBs with 2/3~^ = 0(1) table lookups. □ 



4 Exact Maximum Weight Matching 

At a high level our exact mwm algorithm is similar to our (1 — e)-MWM algorithm. It consists of log TV + 1 
scales, where, in the ith scale, the magnitude of dual adjustments and the violation of domination/tightness 
is bounded in terms of Si, which decreases geometrically with i. However, beyond this similarity the two 
algorithms are quite different. Our exact MWM algorithm only works on bipartite graphs; we assume for 
simplicity that the graph consists of exactly n left vertices and n right vertices. 

Wc redefine Sq = 2L'°s(^/^)J and let Si = Sq/2^ and Wi{e) — Si\_w{e) / Si\ be the granularity and weight 
function of the ith scale, where i e [0, L] and L — [log iV] The algorithm maintains a matching M and 
duals y satisfying Property [3j (As the graph is bipartite there is no need for blossoms or their duals z.) 
Whereas Property [2] allows domination to be violated by Si but enforces tightness of matched edges (of type 
z). Property [3] enforces domination but lets tightness be violated by up to 3^^. 

Property 3. Let i g [0, L] be the scale, M be the current matching, and y : V M>o be the vertex duals, 
where y{e) ~ y{u) + y(v) for edge e = (u, v). 

1. Granularity: y{u) is a nonnegative multiple of Si. 

^^Note that if Af < y/n wc do not require a different weight function at each scale since Wi = w for all i > 0. 
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2. Domination: y{e) > Wi{e) for all e £ E. 

3. Near Tightness: Let e € M be a matched edge. In scale 0, y{e) < u'o(e) + Sq. Throughout scale i we 
have y{e) < Wi{e) + 3(5^ and at the end of scale i we have y{e) < Wi{e) + 5i. 

4. Free Vertex Duals: In scale « = 0, right free vertices have zero y-values and left free vertices have equal 
and minimal y- values among left vertices. At the end of scale and throughout scales 1, . . . , L, all free 
vertices have zero y-values. 

Lemma 10. Let M he the matching at the end of scale i and M* be the MWM. Then w{M) > w{M*) — 2ndi, 
and when i — L, w{M) > ■w{M*) — nS^ > w{M*) - 

Proof. By Property [3] we have 

w{M) > Wi{M) > ^ y{e) — ndi near tightness 

= — nSi free vertex duals 

> Wi{M*) — nSi domination, non-negativity of y 

> w{M*) - 2nS, defn. of Wi, \M* \ < n 

Note that the last inequality is weak whenever 6i < 1, since in this case w = Wi. Specifically, when i — L we 
have Sl = 2Li°s(^/^)J-ri°g^l < l/^, implying w{M) > w{M*) - uSl > w{M*) - □ 

As in our (1 — e)-MWM algorithm we restrict our attention to augmentations on eligible edges. However, 
our definition of eligibility depends on the context. For integers < a < 6, the eligibility graph G[a,b] at 
scale i consists of all edges e such that 

— e ^ M and y{e) = Wi{e), or 

— e G M and Wi{e) + aSi < y{e) < Wi{e) + b5i. 

The algorithm consists of three phases, each with a distinct goal. Phase I, Phase III, and each scale of 
Phase II will require 0{m^/n) time, for a total of 0{my/n\ogN) time. 

Phase I The phase operates only at scale 0. It is a simplified execution of the Gabow-Tarjan [46| algo- 
rithm, stopping not when M is perfect but when y-values of free vertices are zero. In this phase an 
augmentation is an augmenting path in G[l, 1] whose endpoints are free. 



Phase II The phase operates at scales ? = 1, . . . , L. At the beginning of the scale M C G[Q, 3]. The goal is 
to eliminate Af-edges that violate near-tightness by 25i or 35^, that is, the scale ends when M C G[0, 1]. 
In this phase an augmentation is either an augmenting cycle in G[l, 3] or an augmenting path in Gfl, 3] 
whose endpoints have zero y-values. Note that the ends of augmenting paths can be either free vertices 
or matched edges. 



Phase III The phase operates only at scale L. The last scale of Phase II leaves M C G[0, 1]. In this phase 
an augmentation is either an augmenting cycle or an augmenting path whose ends have zero y-values, 
that, in addition, contains at least one non-tight edge. That is, the augmenting path/cycle must exist 
in G[0, 1] but not G[0, 0], which implies that augmenting along such a path increases the weight of the 



matching. By Lemma 10 at the end of Phase II w{M) > w{M*) — ^fn. We guarantee that w(M) is a 
nondecreasing function of time, so, by the integrality of edge weights, k;(M) can be improved at most 
Jn times in Phase III. 
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The fohowing notation will be used liberally. Let X C V{G) be a vertex set, be a subgraph of G, and 
M be an arbitrary matching. Define Vodd{X, H) (and Vcvcn{X, H)) to be the set of vertices reachable from 
X in by an odd-length (and even-length) alternating path starting with an unmatched edge. The directed 
graph H is obtained by orienting e e E{H) from left to right if e ^ A/ and from right to left if e G M . In our 
algorithm H is always chosen to be G[a, h] for some parameters a, h. Note that Vodd(-''', H), V"cvcn(-'^, H), and 
H are defined with respect to a matching known from context. It is clear that VoddiX^H) and Vcvci\{X^H) 
can be computed in linear time, for example, with depth first search (DFS). 

4.1 Phase I 

In Phase I we operate on G[l, 1]. We begin with an empty matching M — % and let y{v) = for right 
vertices and y{v) = So[N/So\ for left vertices. This clearly satisfies Property |3}[2]) (domination) since wo{e) = 
So [w{e)/6o\ and the maximum edge weight is N. The Phase I algorithm oscillates between augmenting along 
a maximal set of eligible augmenting paths and performing dual adjustments. Note that if we augment along 
an eligible augmenting path, all edges of the path become ineligible. See Algorithm [l] for the details. 



ALGORITHM 1: Phase I 

Initialization: 

M ^ 

, , I (5o I A^/(5o I if w is a left vertex 
y(v) ^ < 

10 if u is a right vertex 

repeat 

Augmentation: 

Find a maximal set of augmenting paths in G[l, 1] and set M ^ M ® Upe* ^■ 

Dual Adjustment: 

Let F be the set of left free vertices 

(yiv)-6o if «e Kven(f,G[l,l]) 
y{v) ^ I y{v) +So ifve V,dd{F, G[l, 1]) 
[y(t;) otherwise 
until y-values of left free vertices are zero or all left vertices are matched 



After an augmentation step, there cannot be any augmenting paths in G[l,l], which implies no free 
vertex is in Vodd(^', G[l, 1]). Property |3}|4| is maintained for left free vertices since their y-values are reduced 
in lockstep in every dual adjustment step. It is also preserved for right free vertices since they are not in 
K)dd(-F, G[l, 1]) and therefore never have their y-values adjusted. Property [3j[2]) (domination) is maintained 
since no ehgible edge can have one endpoint in Veven(^7 G[l, 1]) without the other being in Vodd{F, G[l, 1]). 
Property |3]|3| (near tightness) is maintained since for any e S M, y{e) is unchanged if e is eligible, and, if e 
is ineligible (that is, y{e) = wo{e)), y{e) may only be incremented by 6q. 

The number of augmentation/dual adjustment steps is bounded by the number of dual adjustments, that 
is, [S^lN / 5o\) / So < iV/2Li°g(^/v^)J < 2y/n. Thus, Phase I takes 0{my/n) time. 

4.2 Phase II 

At the beginning of scale i e [1, L] in Phase II we set y{u) -s— y{u) + Si for each left vertex u and leave the 
y-values of right vertices unchanged. Since Wi{e) < Wi^i{e) + (5^, this preserves Property |3][2| (domination). 
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Property [sf^S) (near tightness) is also maintained. 



< Wi_i(e) + (5i_i + 5i by Property |3j3) at the end of scale i — 1 

< Wi{e) + 36i since Si^i = 2Si and Wi_i{e) < mi(e) 



However, Property |3j|4|) may be violated since the ?/- values of left free vertices are 6i, not zero. Hence, we 
will run one iteration of Phase I's augmentation and dual adjustment steps on G[l, 3]. These steps preserve 
domination and near tightness and bring left free vertices' y- values down to zero, restoring Property |3j|4| . 
This procedure takes 0{m) time but is executed just once for each of the log scales. 

Next, we will repeatedly perform Phase H Augmentation and Phase H Dual Adjustment steps (see 



Sections 4.3 andK4| on G[l, 3] until M n G[2, 3] = 0, that is y{e) < Wi{e) + di for all e€ M. At this point 



scale i ends and scale i + 1 begins. 



4.3 Phase II Augmentation 

The goal of an augmentation step in Phase II is simply to eliminate all augmenting paths and cycles from 
G[l,3]. We do this in two steps, first eliminating augmenting cycles then paths. Notice that in contrast to 
Phase I, augmenting paths may start or end with matched edges. 

In the first stage of augmentation we will find a maximal set of vertex-disjoint augmenting cycles C using 
DPS. Observe that directed cycles in G[l,3] correspond to augmenting cycles in G[l,3]; this simplifies the 
DPS algorithm since we do not need to distinguish matched and unmatched edges. At all times the DPS 
stack S forms an alternating path. If there is a back edge from the top of the stack to another vertex on this 
stack, that is, S — {. . . ,v, . . . ,u) and (u, v) e E{G[1, 3]), this back edge closes an augmenting cycle, which 
can be added to C. A vertex is marked when it is popped off the DPS stack, either because the vertex joins 
an augmenting cycle in C or if it is found not to be contained in any augmenting cycle. Algorithm [2] gives 
the details for Cycle-Search. 

Lemma 11. The algorithm Cycle-Search finds a maximal set of vertex-disjoint augmenting cycles C. More- 
over, if we augment along every cycle in C , then the graph G[l, 3] contains no more augmenting cycles. 

Proof. Supposing C is not maximal, let G = (vq,vi, . . .Vk-i,vo) be any cycle vertex-disjoint from all cycles 
in C, where Vq is the first vertex of G pushed onto the stack. Let t be the largest index such that Vt is pushed 
onto the stack before Vq is popped off. It follows that vt+inwdk is unmarked and therefore appears in Vy^. 
If f = fc — 1 then Vq € Vvf,_-^ and the search will discover an augmenting cycle containing ut; if t < fc — 1 then 
the search will push Vt+i onto the stack, contradicting the maximality of t. 

If there exists an eligible cycle G after augmentation, then this cycle must share a vertex v with some 
cycle G' e C due to the maximality of C. However, since G' contains v's mate (both before and after 
augmentation), G and G' must intersect at an edge, which contradicts the fact that all edges in G' become 
ineligible after augmentation. □ 

In the second stage of augmentation wc will eliminate all the augmenting paths in G[l, 3]. This is done by 
finding a maximal set of vertex-disjoint maximal augmenting paths, which are those not properly contained 
in another augmenting path. (Recall that in Phase II, augmenting paths can end at matched edges so long 
as the endpoints of the path have zero y-values. Any augmenting path with free endpoints is necessarily 
maximal, but one ending in matched edges may not be maximal.) See Figure [s] for an illustration. 

Consider the graph G[l,3]. It must be a directed acyclic graph, since, by Lemma [n] G[l,3] does not 
contain augmenting cycles. Call a vertex with zero y-value a starting vertex if it is left and free or right 
and matched, and an ending vertex if it is left and matched or right and free. Let A and B be the set of 
starting and ending vertices. It follows that any augmenting path in G[l,3] corresponds to a directed path 
from an A- vertex to a _B-vertex in G[l,3]. With this observation in hand we can find a maximal set V of 
maximal augmenting paths using DPS. We initiate the search on each A-vertex u in topological order. When 
an augmenting path to a B- vertex, say v, is first discovered we cannot commit (u, . . . ,v) to V immediately. 
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ALGORITHM 2: Cycle- Search: returns a maximal set of eligible augmenting cycles 
All vertices are initially unmarked 

C % {set of augmenting cycles discovered so far} 

5 -s- {DFS stack} 

while there are still unmarked vertices do 

if S is empty then 

Push any unmarked Uq &V onto S. 

end if 

u = the top of S 

Ki = {w I (m, v) e _E(G[1, 3]) and v is unmarked} 
if Vu is empty then 

Mark u and pop u off S 
else 

Choose any v &Vu 

if V appears in S, that is, if 5 = (. . . , . . . , m) then 

C C U {(w, . . . ,u, v)} {add augmenting cycle to C} 

Mark u, . . . , m and pop (v, . . . ,u) off 5 
else 

Push w on 5 
end if 
end if 
end while 
return C 



Rather, we continue to look for an even longer augmenting path, and keep {u, . . . ,v) only after all outgoing 
edges from v have been exhausted. Algorithm [3] gives the Path-Search procedure. 

Lemma 12. After augmenting along every path in V, the graph G[l,3] contains no augmenting paths. 

Proof. Suppose that there exists an augmenting path Q after the augmentation. Then, by the maximality 
of V, there must be some augmenting path in V sharing vertices with Q. There can be two cases. 

Case 1 There exists a P G V and a v G P D Q that is not an endpoint of P. Since P contains v and its 
mate before and after augmentation, P and Q must share an edge, which is impossible since all edges 
in P become ineligible after augmentation. 

Case 2 Let P be the first path added to V that intersects Q. Since we are not in Case 1, P and Q intersect 
at a common endpoint, say x. (Note that the only way this is possible is if x is matched before 
augmentation and free afterward. If x were a free endpoint of P before augmentation it would be 
incident to a matched ineligible edge after augmentation and could therefore not be an endpoint oi Q.) 
Let xp and xq be the other endpoints of P and Q. If xp — xq then G[l, 3] contained an augmenting 
cycle, contradicting Lemma If G A is a starting vertex, consider the moment when the stack 
contained only xp. At this time PQ is an augmenting path of unmarked vertices, so Path-Search would 
not place the non-maximal augmentation P in V . On the other hand, if a; is a starting vertex, it must 
be a right matched vertex that becomes free after augmentation, which implies that xq must also be 
a starting vertex. Since our search explores A-vertices in topological order, the search from xq would 
have preceded the search from x. Thus, the first augmentation of V intersecting Q must contain xg, 
contradicting the fact that P does not contain xq. 

□ 
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Figure 3: An illustration of starting vertices and maximal augmenting paths in Gil, 3]. The plain edges 
denote unmatched edges; the shaded edges are matched. The haloed vertices have zero y-values. The 
set of starting vertices is {ui,vi,V2}- The path P — (w2, U3, W3, U4, W4) is an augmenting path but not a 
maximal augmenting path, since it can be extended to a longer one. For example, (wi, W2, W3, W3, 1*4, W4) and 
{vi, U2,V2,U3,V3,U4, V4) are augmenting paths containing P. 



4.4 Phase II Dual Adjustment 

Recall that scale i ends when M C G[0, 1]. Let B ^ Mn G[2, 3] be the bad edges and let f : E ^ {0, 1, 2} 
measure the badness, defined as follows. 



/(e) 



' for e ^ B 



Let /(M) = J2eeB /(e) be the total badness of M. The goal of dual adjustment is to eliminate B, or 
equivalently, to reduce f{M) to 0. We will show how to reduce /(M) by roughly f(M) in linear time. 
A B' <Z B \s called a chain if there is an alternating path in G[l,3] containing B' and an anti-chain if 



no alternating path in G[l,3] contains two edges 61,62 G B' . Lemma 13 basically follows from Dilworth's 
lemma |22| . 

Lemma 13. For any t > \, there exists B' Q B such that either B' is a chain with J{B') > \t~\ or B' is an 
anti-chain with \B'\ > \ f{M)/2t~\ . Moreover, such a B' can be found in linear time. 

Proof. Let S be the set of vertices with zero in-degree in the acyclic graph G[l, 3]. In linear time we compute 
distances from S using —/as the length function. Let d{v) be the distance to v. Suppose there is a u with 
^ — r^l and let P be a shortest path to v. It follows that B' — P D B is a. chain with f{B') > [t] . If this 
is not the case then, for every {u,v) G B (where w is a left vertex), d{v) G [—{\t'] — 1), —1]. Since /(e) < 2 for 
6 G -B, we must have at least [|i?|/([<] — 1)] > [/(M)/2i] such v with a common distance, say —k. It follows 
that B' = {{u,v) G -B I is a left vertex and d{v) — —k} is an anti-chain, for if {ui,vi), (^2,^2) G B' were 
on an alternating path in G[l, 3], the distance to V2 would be strictly smaller than the distance to vi. □ 

Below we show that if B' is a chain we can decrease the total badness by f{B') in linear time. On the 
other hand, if B' is an anti-chain, then we can decrease the total badness by \B'\/2, also in linear time. 



4.4.1 Phase II Dual Adjustment: Antichain Case 

When performing dual adjustments we must be careful to maintain Property [sj which states that y{u) must 
always be nonnegative, and must be zero if u is free. This motivates the definition of a dual adjustable vertex. 

Definition 4. A vertex u is said to be dual adjustable if every vertex in V^dd(^*j G[l, 3]) is matched and 
every vertex v G Vi3von(u, G[l, 3]) has y{v) > 0. 
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ALGORITHM 3: Path-Search: returns a maximal set of maximal eligible augmenting paths 
All vertices are initially unmarked 

A -(^ {u\ y{u) = and u is left and free or right and matched} {starting vertices} 

B ^ {u I y(u) — and u is right and free or left and matched} {ending vertices} 

T' -s— {set of augmenting cycles discovered so far} 

5^0 {DFS stack} 

while there are still unmarked A-vertices do 

if S is empty then 

Push the first (topologically) unmarked wq € A onto S. 

end if 

u = the top of S 

yu~{v\ {u, v) £ E{G[1, 3]) and v is unmarked} 
if Vu is empty then 
if u E B then 

V ^ V U {S} {S forms a maximal augmenting path} 

Mark all vertices in S and set S ^ {) {Pop all vertices off the stack} 

else 

Mark u and pop u off 5 
end if 
else 

Push an arbitrary v E Vu onto S 
end if 
end while 
return V 



Lemma 14. For every e ~ {u, v) e B, either u is adjustable or v is adjustable or both. Furthermore, all 
adjustable vertices can be found in 0{m) time. 

Proof. First, suppose that for some e = (m, v) G B, both u and v are not adjustable. There must exist 
vertices w and x having zero y- values where (w, . . . ,u,v, . . . ,x) is an augmenting path in G[l, 3]. However, 
this contradicts Lemma |l2[ which states that there are no augmentations in G[l,3] after an augmentation 
step. Let V — {v \ V is free or {v, v') € M and y{v') — 0}. By definition a vertex is not dual adjustable if 
and only if it lies in Vodd(^j 3]), which can be computed in linear time. □ 

Let B' C B he an anti-chain. The procedure Antichain-Adjust{B') selects a set of dual adjustable vertices 
X incident to B' and on a common side (left or right), then does a dual adjustment starting at X. Since, 



by Lemma 14 for any (m, v) € B' either u is adjustable or v is adjustable or both, we can guarantee that 
1^1 > \B'\/2. See Figure|4]for an example. 

Lemma 15. The dual adjustment starting at X will not break Property 3- 1\S\4 )■ Furthermore, it makes 
PropertylMSy tighter by decreasing f{M) by \X\. 



Proof. Since X consists of adjustable vertices, every vertex v S VcvoniX, G[l, 3]) must have y{v) > 0, implying 
y{v) will be non-negative after being decremented by Si. Thus, Propert y [3p| ) is maintained. Furthermore, 
Vodd{X, G[l, 3]) cannot contain a free vertex, which implies that Property |3[|4[) is preserved. Since all vertices 
in X are on the same left/right side, y(e) can change by at most Si. PropertyT3l[2]) (domination) is maintained 
since no tight edge (with y(e) = ^^(e)) can can have one endpoint in Vcvcn(A, G[l, 3]) without the other being 
in yodd(-''^, 3]). The algorithm does not increase f{M) since, for any e S MnG[l,3], if one endpoint of e 
appears in Vodd(A, G[l, 3]), the other must appear in V'cvcn(A, G[l, 3]). Furthermore, the algorithm decreases 
f{M) by at least \X\ > \B'\/2 since f{u,v) is decremented for each {u,v) e B' with v £ X. To see this, 
note that y{v) is decremented by Si and, since G[l,3] is acyclic and B' is an antichain, u cannot appear in 
ViDdd(^, G[l, 3]) and therefore cannot have its y-value adjusted. □ 
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ALGORITHM 4: Antichain-Ad]ust{B') 



Find adjustable vertices: 

V -f^ {v\ V is free or (w, v') g M and y{v') — 0}. 

Mark vertices in \ V;3dd(V^, 3]) as adjustable. 

Xl {u I (u, w) e -B' and u is a left adjustable vertex}, 

Xr ^ {u I (u, w) G B' and m is a right adjustable vertex}. 

If \Xr\ > \Xl\, then let X ^ Xr; otherwise let X ^ Xl- 



Perform dual adjustments, starting at X: 

'yiv)^S, if yeve„(X,G[l,3]) 
y{v)^L{v) + 6, if«GFodd(X,G[l,3]) 
y{v) otherwise 




(a) 




(b) 



Figure 4: (a) The haloed vertices have zero y-values. The circled matched edges form an antichain B' . The 
shaded vertices in V(B') are dual adjustable, (b) The black vertices are in X and vertices marked with 'e' 
and 'o' are in Vcvcn{X,G[l,3]) and T4dd(-'^, 3]), respectively. 



Therefore, by doing the dual adjustment starting at X, we can decrease /(M) by at least \B'\/2. 
4.4.2 Phase II Dual Adjustment: Chain Case 

In the chain case we are given a chain B' C B and a minimal alternating path P containing B' , that is, it 
starts and ends with B'-edges. Setting M ^ M © P immediately reduces f{M) by f{B') since P'-edges 
are replaced by tight edges, which contribute nothing to f{M). However, the endpoints of P, say u and w, 
are now free while possibly having positive y- values, which violates Property [3j[4]) . Our goal is to restore 
Property [Sj either by finding augmenting paths that rematch u and v or by reducing their j/- values to zero. 
In this section an augmenting path is one in the eligibility graph G[0, 3] (not G[l, 3]) that either connects u 
and V or connects an a; € {u, v} to a vertex with zero y-value. A notable degenerate case is when y{x) = 0, 
in which case the empty path is an augmenting path from x to x. We begin by performing dual adjustments 
(as in a Hungarian search) until an augmenting path P„ in G[0, 3] containing u emerges. We do not augment 
along P[j immediately but perform a second search for an augmenting path P^, containing v. If P« and P„ 
do not intersect we let Q = PuU Py. On the other hand, if Pu and P^, do intersect then there must be an 
augmenting path P„„ between u and v in G[0, 3]; we let Q = P,„,. In either case we augment along Q, setting 
M ^ M (B Q- Figure [5] illustrates the case where Q — Puv 

The search from x € {u, v} works as follows. If there exists an augmenting path P^ in G[0, 3] starting at 
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X, then return P^. (This will be the empty path if y{x) = 0.) Otherwise update y- values as follows: 



ly{z)-S, ifzG V;vc„(2:,G[0,3]) 
viz) ^ I y{z) + 5, if z e V;dd(a:, G[0, 3]) 
I y(z) otherwise 

and continue to perform these dual adjustments until an augmenting path starting from x emerges. As 
in standard Hungarian search, this process is reducible to single source shortest paths on a non-negatively 
weighted directed graph. The graph is either G or its transpose, depending on the left/right side of x. The 
weight function w is zero on M-edges and equal to the gap between y-value and weight on non-M edges: 
w{e) = y{e) — Wi{e). By Property [3j[2]) (domination), w is non-negative. The value h{z) is the sum of dual 
adjustments until an augmenting path emerges from x to z. If z is free and on the opposite side of x then 
this is simply the distance from x to z. However, if z is on the same side as x (and therefore matched), 
we need another y{v)/5i dual adjustment steps to reduce its y-value to zero. The pseudocode for Search{x) 
appears below. 



ALGORITHM 5: Search{x) 



G< 
'w{e) <— 



Initialize the weighted graph: 

G if x is a left vertex 

G^ if x is a right vertex (reverse the orientation) 

y{e) - Wi{e) if e ^ M 
^0 ifeeA/ 
Compute distances. For each z ^ V: 

d{z) distance from a; to z in G with respect to w, or oo if z is unreachable from x. 

{d{z) if z is free and not on the same side as x 

d{z) + y{z) if z is on the same side as x 
oo otherwise 
2^inin arg min^g y{/i(z)} 

A ^ /l(Zniin) 

Perform dual adjustments. For each z G V: 

^ ^ I y{z) — max{0, A — d{z)} if z is on the same side as x 

1 y{z) + max{0, A — d{z)} if z is not on the same side as x 
Pj; a shortest path from x to Zmin- (Note: 'w{Px) = after dual adjustments above.) 
return 



Each execution of Search{x) clearly takes linear time, except for the computation of the distance function 
d. We use Dijkstra's algorithm [3T], implementing the priority queue as an array of buckets [20]. Since 
'w{e)/5i is an integer and we are only interested in distances at most A (see the pseudocode of Search), 



a ( A/ (5i) -length array suffices. Lemma 16 implies that A = 0(n6i), which gives a total running time of 
0(to -|- n). 

Lemma 16. Augmenting along P then Q does not decrease the weight of the matching, that is, w{{M © 
P) (B Q) > w{M). Furthermore, A„ + A^, < 3nSi, where is the sum of dual adjustments performed by 
Search(x). 

Proof. Call Ml = M (B P and M2 = Mi © Q the matchings after augmenting along P and then Q and let 
w be the weight function w;(e) = y{e) — Wi{e). (Notice that w differs from w on the matched edges.) For a 
quantity q denote its value before Search{u) and Search{v) by qoid and after both searches by qnew After 
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y{u) = 10(5. 



O- 



3(5,: 



(a) 



O O 



O V 



25, 

y{v) = m, 



U O 





y{u) = 6Si 



Pu 



o- 



(b) 



-O 



85, 



O 



O V 



26, 

y{v) = 10<5, 



UQ 







=0 



o 



65,; 





y{v) = 6(5,: 



O V 



u O 



y{u) = 6(5, 



O 



o- 



-O 



=0 V 



y{v) = 6(5, 



(c) 



(d) 



Figure 5: An example illustrating procedures for the chain case. Edges are shown with their new weight w. 
The haloed vertices are free vertices with zero y-values. (a) After augmenting along P, u and v become free 
while having positive ?/- values, (b) search{u) adjusts duals by A„ = A5i and finds an augmenting path P„. 
(c) search{v) adjusts duals by A„ = A5i and finds P^. (d) Augmentation along Q. This is the case where 
there exists an augmenting path Q between u and v in G[Q, 3], which happens to be P^, in the example. 



the two searches, we must have: 



,(Q\Mi)= Vnewie) 
eeQ\Mi 

= Vnewiu) + Vnewiv) + ^ Vnewis) 
eeMinQ 

— Unew 

{Ml n Q) 



tightness on unmatched edges 



defn. of Wn 



(6) 



Line (|6| follows since, aside from u and v, V{Q \ Mi) and V{Mi n Q) differ only on vertices with zerc 
y- values. (These are the other endpoints of P^ and Py when Q = P^U Py .) Therefore, 



w,{M2) = w.,{Mi) + w,(Q \ Ml) - w,{Mi n Q) 

= Wi{Mi) + (Ml n Q) 

A similar proof shows that before the two searches, we have 



(7) 



w,(M) = w,(Mi) + yoid{u) + Voidiv) - Woid{M n P) (8) 
The total amount of dual adjustment in both searches is at most the distance from u to v, so Au + A„ < 



Woid{M n P) < 3ndi. Moreover: 

Wi{M2) > Wi{Mi) + y„ew{u) + Vnewiv) 

= w,(Mi) + yoid{u) + Voidiv) - Ay- Ay 
> u;,(Mi) + yoid{u) + yoid{v) - Woid{M n P) 



by and Wnew{Mi n Q) > 
by (§ 



□ 
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Lemma 17. At most 0{y/n) rounds of augmentation and dual adjustment are required to reduce f{M) to 
0. 



Proof. We apply Lemma 13 with t = Vb/2, where b = f{M). It follows that in linear time we can either 
obtain an anti-chain B' of size at least [v^] and reduce f{M) by [•\/6/2], or obtain a chain B' such that 
f{B') > \Vb/2'\ and reduce /(M) by f{B'). In either case we can reduce /(M) by \Vb/2\. The number 
of rounds is at most T{b), where T(0) = and T{b) = T{b - \Vb/2']) + 1 for 6 > 0. A proof by induction 



shows that T(6) < 4\/6, which is at most 4%/2n since \M\ < n and /(e) < 2 for e G M. □ 

Phase II concludes at the end of scale L = [log N~\ , when Mn G[2, 3] — 0. That is, when Phase III begins 
y{e) — w{e) or w{e) + Sl for each e G Af. (Note that wl = w since (5l < 1.) 



4.5 Phase III 

Like a single scale of Phase II, Phase III alternately executes augmentation and dual adjustment steps. 
Certain complications arise since the goal of Phase III is to eliminate all non-tight M-edges whereas the 
goal of Phase II was to eliminate edges in M n G[2,3]. In order to understand the ramifications of this 
slight shift, we should review the interplay between augmentation and dual adjustment in Phase II. Note 
that Phase II augmentation steps improved the weight of M, but this was an incidental benefit. The real 
purpose of augmentation was to eliminate any augmentations in G[l,3] (thereby making G[l,3] acyclic), 
which then let us reduce f{M) by ^ f(M)/2 with a chain/ ant ichain dual adjustment. The efficiency of 
the augmentation step stemmed from the fact that matched and unmatched edges had different eligibility 
criteria. Thus, augmenting along a maximal set of augmentations in G[l,3] destroyed all augmentations in 
G[l,3], that is, tight M-edges could and should be ignored. 

In Phase III we cannot afford to exclude tight M-edges from the eligibility graph. This raises two 
concerns. First, augmenting along a maximal set of augmenting paths/cycles in G[0, 1] does not destroy all 
augmentations in G[0, 1]. (Tight M-edges remain eligible after augmentation and may therefore be contained 
in another augmentation.) Second, G[0, 1] may contain cycles of tight edges, that is, augmentations that 
do not improve the weight of M , so eliminating all weif^/if- augmenting paths/cycles does not guarantee that 
G[0, 1] is acyclic. These concerns motivate us to redefine augmentation. In Phase III an eligible augmentation 
is either an alternating cycle or alternating path whose endpoints have zero y-values, that, in addition, is 
contained in G[0, 1] but not G[0, 0]. That is, it cannot consist solely of tight edges. 



4.5.1 Phase III Augmentation 

In a Phase III Augmentation step we repeatedly augment along eligible augmentations in G[0, 1] until no 



such eligible augmentation exists. Lemma [18] lets us upper bound the aggregate time for all Phase III 
Augmentation steps. 

Lemma 18. An eligible augmentation P in G[0, 1] can be found in 0{m) time, if one exists, and w{M®P) > 
w{M). Consequently, there can be at most ^Jn augmentations in Phase III. 

Proof. To find an augmentation we begin by computing the strongly connected components of G[0, 1] in 
linear time. If the endpoints of an e G M n G[l, 1] are in the same strongly connected component then that 
edge is contained in an augmenting cycle. If there are no augmenting cycles, use the linear time algorithm 



described in Lemma 14 to determine whether v is dual adjustable for all v £ V. If both endpoints of some 
e G M n G[l, 1] are not dual adjustable, then there must be an augmenting path containing e. 

Let P be an augmenting path or cycle. We must have X]«iEy(pnAf) i'(^) ~ J2vev{P\M) vi'")^ since the sums 
differ only on vertices with zero y- values. Thus w(PnM) < J2vev(PnM) V^"^^ ~ '^vev(P\M) vi"^) ~ 'w{P\M), 
where the first inequality follows from the fact that POM contains a non-tight edge and the last equality 
from the tightness of unmatched edges. Since all weights are integers, the weight of the matching is increased 
by at least one. 
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Recall that 5l = 2Li°g^/V^J-ri°g^'l < l/^. By Lemma[lo) we have w{M) > w{M*) - uSl > w{M*) - 
y/n at the end of Phase II, where A/* is the mwm. As in Phase II, Phase III dual adjustments will not 



reduce the weight of M. (See Lemma 16 ) Since each augmentation increases the weight of M there can be 
at most ^/n augmentations in Phase III. □ 



Lemma 



18 



implies that the time for £ Phase III Augmentation steps is 0{£m + m^/n): linear time per 



step plus linear time per augmentation discovered. 



4.5.2 Phase III Dual Adjustment 



Since the goal of Phase III is to eliminate edges in M n G[l, 1] (as opposed to M n G[2, 3]) we must redefine 
the badness of an edge accordingly. Let B — M O G[l, 1] be the bad edges and let f : E ^ {0, 1} measure 
badness, where /(e) = {y{e) — w{e))/5L if e € B and zero if e ^ i3. We define B' Q B io he & chain or 
antichain exactly as in Phase II. The difference is that G[0, 1] is not necessarily acyclic so finding a chain or 
antichain requires one extra step. 

Lemma 19. For any t > 1, there exists a B' Q B such that B' is a chain with f{B') > \t~\ or B' is an 
antichain with \B'\ = f{B') > \f{M)/t~\. Moreover, B' can he found in linear time. 

Proof. Let G* be the graph obtained from G[0, 1] by contracting all strongly connected components in 
G[0, 0]. Since G[0, 1] contains no eligible augmenting cycles after a Phase III Augmentation step, all i3-edges 
straddle different strongly c onn ected components and therefore remain in G*. By definition G* is acyclic. 



The argument from Lemma 13 shows that G* contains a chain B' with f{B') > \t~\ or an antichain with 



B'\ — \f{AI)/t'] , and that such a B' can be found in linear timej^ 



□ 



We can apply the chain and antichain elimination procedures from Phase II without compromising 
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15 



16 



and 
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remain valid if we substitute G[0, 1] for G[l, 3]. Let t 
Lemmas 



6/2, 



14 



and 19 imply that we can reduce f{M) by 



correctness since Lemmas 

where b = f{M) is the current total badness. 

\t~\ > \fbj2 in the chain case or [6/2f] > \/b/2 in the antichain case. The number of augmentation and dual 
adjustment steps in Phase III is then T{b) where T{b) = T{b - \\/b/2']) + 1 and T(0) = 0. By induction 
T{b) < 2\/26, which is at most 2^/2n since /(M) < n. Thus, the total time spent on dual adjustment in 
Phase III is still 0{m^n). 



4.6 Maximum Weight Perfect Matching 

Recall from Section [l] that the maximum weight perfect matching problem (mwpm) is reducible to mwm. 
One simply adds nN to the weight of every edge; a mwm in the new graph is necessarily a mwpm in the 
original. Thus, our mwm algorithm solves the mwpm problem in 0(r7iY^log(niV)) time, where the number 
of scales is log((n+l)A^). However, we can circumvent this roundabout reduction and solve mwpm directly, in 
log(Y^A^) scales, that is, a factor 2 improvement for small N n. We substitute Property [4] for Propertyjs] 

Property 4. Redefine Sq = 2L'°s^J, L — \log{^y nN)'] . In each scale i € [0, -Z^], we maintain a perfect 
matching M satisfying the following properties. 

1. Granularity: y{u) is a multiple of 5i for all u V. 

2. Domination: y(e) > Wi{e) for all e d E. 

3. Near Tightness: For any e £ M, y[e) < Wi{e) + 3Si throughout scale i and y{e) < Wi{e) + Si at the 
end of scale i. 

i^The bound on \B'\ is [/(M)/i] rather than [/(M)/2t] since the range of / is {0, 1} rather than {0, 1, 2}. 
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yiu) 



In Phase I we find any perfect matching M in 0{m^/n) time [551 [531 [73] and initialize values to satisfy 
Property [4] 

(5o if u is a left vertex 
if u is a right vertex 

Since j/(e) < Wi{e) + for e e Af , Property |4] lets us end scale 0. 

As before, Phase II operates at scales i £ [i,L] and Phase III operates at scale L with the following 
simplifications: 

1. To begin scale i we simply increment y{u) by Si for each left vertex u. There is no need for an initial 



augmentation/dual adjustment (as in Section 4.2 1 since there are no free vertices 



2. As there are no free vertices, an augmentation is always an augmenting cycle. In the Phase II and 
Phase III augmentation steps we use only Cycle-Search, not Path-Search. 

3. We do not prohibit negative y-values. Thus, in the antichain case of dual adjustment, both endpoints 
of a -B'-edge are dual-adjustable and we can reduce /(M) by \B'\ rather than just \B'\/2. 

4. In the chain case of dual adjustment, setting M •(— M ® P temporarily frees P's endpoints, say u 
and V. We execute Search{u) but force z„iin = v, that is Search returns a u-to-v path P„. Setting 
M ^ AI (BPu restores the perfection of M. There is no need for Search{u) to calculate /i- values. These 
were introduced to maintain Property ^HTl), namely that y- values are non- negative. 



We change the parameter that determines whether we apply a chain or antichain dual adjustment. 
Choose t = ■y/5/2, where b = f{M). Either we can obtain an anti-chain B' of size at least [■y/6/2] 
and decrease f{M) by [■y/6/2] , or we can obtain a chain B' such that f{B') > \y/b/2] and decrease 
f{B) by f{B'). Thus, the number of rounds is at most T{b) = T{b - \Vb/2\) + 1 < 2\/26. In a Phase 
II scale h < 2n, so Ay^n augmentation/dual adjustment steps are needed. In Phase III 6 < n so 2y/2n 
steps are needed. 

Note that since 6l = 2Li°g^J-ri°g(v^^)l < l/^n, at the end of Phase II we have w{M) > w{M*) - 



tiSl > w(M*) — by Lemma 10 Thus, the matching can be augmented at most ^/n times in Phase 



III and the total time spent in augmentation steps is 0{my/ri). 

The main difference between our mwm and mwpm algorithms is in Phase I. In the mwm algorithm we 
can afford to use a smaller value for 5q since Phase I ends when free vertices have zero y-values, whereas 
Phase I of the mwpm algorithm ends only when we have a perfect matching. 



References 

[1] p. K. Agarwal and K. R. Varadarajan. A near-linear constant-factor approximation for Euclidean 
bipartite matching? In Proceedings 20th ACM Symposium on Computational Geometry, pages 247- 
252, 2004. 

[2] H. Alt, N. Blum, K. Mehlhorn, and M. Paul. Computing a maximum cardinality matching in a 
bipartite graph in time 0{n^-^ ^Jm/ \ogn). Info. Proc. Lett, 37(4):237-240, 1991. 

[3] T. Anderson, S. Owicki, J. Saxe, and C. Thacker. High speed switch scheduling for local area networks. 
ACM Trans. Comput. Syst., ll(4):319-352, 1993. 

[4] A. Andersson, T. Hagerup, S. Nilsson, and R. Raman. Sorting in linear time? J. Comput. Syst. Sci., 
57(l):74-93, 1998. 

[5] D. Avis. Two greedy heuristics for the weighted matching problems. In Proceedings 9th Southeast 
Conference on Combinatorics, Craph Theory, and Computing ( Congr. Numer. XXI), pages 65-76, 
1978. 



29 



[6] D. Avis. A survey of heuristics for the weighted matching problem. Networks, 13:475-493, 1983. 

[7] M. L. Balinski. Labelling to obtain a maximum matching. In R. C. Bose and T. A. Downing, editors. 
Combinatorial Mathematics and its Applications, pages 585-602. University of North Carolina Press, 
1969. 

[8] M. L. Balinski and R. E. Gomory. A primal method for the assignment and transportation problems. 
Management Science, 10(3):578-593, 1964. 

[9] D. .P. Bertsekas. A new algorithm for the assignment problem. Mathematical Programming, 21:152- 
171, 1981. 

[10] G. Birkhoff. Tres observaciones sobre el elgebra lineal. Universidad Nacional de Tucuman, Revista A, 
5(1-2):147-151, 1946. 

[11] H. E. Brogden. An approach to the problem of differential prediction. Psychometrika, 11 (3): 139-154, 
1946. 

[12] G. Brown and J. von Neumann. Solutions of games by differential equations. In H. Kuhn and A. Tucker, 
editors, Contributions to the Theory of Games, volume 24 of Annals of Mathematical Studies, pages 
73-79. Princeton University Press, 1950. 

[13] M. Charikar. Similarity estimation techniques from rounding algorithms. In Proceedings 34th ACM 
Symposium on Theory of Computing (STOC), pages 380-388, 2002. 

[14] J. Cheriyan and K. Mehlhorn. Algorithms for dense graphs and networks on the random access 
computer. Algorithmica, 15(6) :521 549, 1996. 

[15] N. Christofi-des. Worst case analysis of a new heuristic for the travelling salesman problem. Technical 
report, Graduate School of Industrial Administration, Carnegie Mellon University, 1976. 

[16] R. Cole, R. Hariharan, M. Lewenstein, and E. Porat. A faster implementation of the Goemans- 
Williamson clustering algorithm. In Proceedings 12th ACM-SIAM Symposium on Discrete Algorithms 

(SODA), pages 17-25, 2001. 

[17] W. H. Cunningham and A. B. Marsh, III. A primal algorithm for optimum matching. Mathematical 

Programming Study, 8:50- 72, 1978. 

[18] G. B. Dantzig. Application of the simplex method to the transportation problem. In T. C. Koopmans, 
editor. Activity Analysis of Production and Allocation, Cowles Commission Monograph 13, pages 359- 
373. John Wiley and Sons, New York, 1951. 

[19] J. F. Desler and S. L. Hakimi. A graph-theoretic approach to a class of integer-programming problems. 

Operations Research, 17(6): 1017-1033, 1969. 

[20] R. B. Dial. Algorithm 360: shortest-path forest with topological ordering. Comm. ACM, 12(11):632- 
633, 1969. 

[21] E. W. Dijkstra. A note on two problems in connexion with graphs. Numerische Mathematik, 1:269-271, 
1959. 

[22] R. P. Dilworth. A decomposition theorem for partially ordered sets. The Annals of Mathematics, 
51(1):161-166, 1950. 

[23] E. A. Dinic. Algorithm for solution of a problem of maximum flow in networks with power estimation. 
Soviet Math. Dokl, 11:1277-1280, 1970. 



30 



[24] E. A. Dinic and M. A. Kronrod. An algorithm for the solution of the assignment problem. Soviet 
Math. Dokl, 10(6):1324-1326, 1969. 

[25] D. Drake and S. Hougardy. Improved linear time approximation algorithms for weighted matchings. 
In Ith International Workshop on Randomization and Approximation Techniques in Computer Science 
(APPROX), LNGS 2764, pages 14-23, 2003. 

[26] D. Drake and S. Hougardy. A simple approximation algorithm for the weighted matching problem. 
Info. Proc. Lett., 85:211-213, 2003. 

[27] R. Duan and S. Pettier Approximating maximum weight matching in ncar-lincar time. In Proceedings 
51st IEEE Symposium on Foundations of Computer Science (FOCS), pages 673-682, 2010. 

[28] I. S. Duff and J. R. Gilbert. Maximum-weighted matching and block pivoting for symmetric indefinite 
systems. In Householder Symposium XV Book of Abstracts, pages 73-75, 2002. 

[29] T. E. Easterfield. A combinatorial algorithm. J. London Math. Soc, 21:219-226, 1946. Republished 
as: An algorithm for the allocation problem. Operations Research 11, 3, pp. 123-129, 1960. 

[30] J. Edmonds. Maximum matching and a polyhedron with 0, 1-vertices. J. Res. Nat. Bur. Standards 
Sect. B, 69B:125-130, 1965. 

[31] J. Edmonds. Paths, trees, and flowers. Canadian Journal of Mathematics, 17:449-467, 1965. 

[32] J. Edmonds and E. L. Johnson. Matching, Euler tours, and the Chinese postman. Mathematical 
Programming, 5:88-124, 1973. 

[33] J. Edmonds and R. M. Karp. Theoretical improvements in algorithmic efficiency for network flow 
problems. J. ACM, 19(2):248-264, 1972. 

[34] T. Feder and R. Motwani. Clique partitions, graph compression and speeding-up algorithms. J. Com- 
put. Syst. Sci., 51(2):261-272, 1995. 

[35] L. R. Ford and D. R. Fulkerson. Flows in Networks. Princeton University Press, 1962. 

[36] M. L. Fredman and R. E. Tarjan. Fibonacci heaps and their uses in improved network optimization 
algorithms. J. ACM, 34(3):596-615, 1987. 

[37] M. L. Fredman and D. E. Willard. Surpassing the information-theoretic bound with fusion trees. 
J. Comput. Syst. Sci., 47(3):424-436, 1993. 

[38] H. N. Gabow. Implementation of algorithms for maximum matching on nonbipartite graphs. Ph.D. 
thesis, Stanford University, 1974. 

[39] H. N. Gabow. Scaling algorithms for network problems. In Proceedings 24th IEEE Symposium on 

Foundations of Com,puter Science (FOCS), pages 248- 257, 1983. 

[40] H. N. Gabow. A scaling algorithm for weighted matching on general graphs. In Proceedings 26th IEEE 
Symposium on Foundations of Computer Science (FOCS), pages 90-100, 1985. 

[41] H. N. Gabow. Scaling algorithms for network problems. J. Comput. Syst. Sci., 31(2):148-168, 1985. 

[42] H. N. Gabow. Data structures for weighted matching and nearest common ancestors with linking. 
In Proceedings First Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 434-443, 
1990. 

[43] H. N. Gabow, Z. Galil, and T. H. Spencer. Efficient implementation of graph algorithms using con- 
traction. J. ACM, 36(3):540-572, 1989. 



31 



[44] H. N. Gabow and S. Pcttic. The dynamic vertex minimum problem and its application to clustering- 
type approximation algorithms. In Proceedings 8th Scandinavian Workshop on Algorithm Theory 
(SWAT), pages 190-199, 2002. 

[45] H. N. Gabow and R. E. Tarjan. Almost-optimum speed-ups for bipartite matching and related prob- 
lems. In Proceedings 20th Annual ACM Symposium on Theory of Computing (STOC), pages 514-527, 
1988. 

[46] H. N. Gabow and R. E. Tarjan. Faster scaling algorithms for network problems. SIAM J. Comput., 

18(5):1013 1036, 1989. 

[47] H. N. Gabow and R. E. Tarjan. Faster scaling algorithms for general graph-matching problems. 
J. ACM, 38(4):815-853, 1991. 

[48] Z. Galil, S. Micali, and H. N. Gabow. An 0{EV logV) algorithm for finding a maximal weighted 
matching in general graphs. SIAM J. Comput, 15(1):120-130, 1986. 

[49] P. Giaccone, E. Leonardi, and D. Shah. On the maximal throughput of networks with finite buffers 
and its application to buffered crossbars. In Proceedings 24th INFOCOM, volume 2, pages 971-980, 
2005. 

[50] A. Gleyzal. An algorithm for solving the transportation problem. J. Res. Nat. Bur. Standards, 
54(4):213-216, 1955. 

[51] M. X. Goemans and D. P. Williamson. A general approximation technique for constrained forest 
problems. SIAM J. Comput, 24(2):296-317, 1995. 

[52] A. V. Goldberg. Scaling algorithms for the shortest paths problem. SIAM J. Comput, 24(3):494-504, 
1995. 

[53] A. V. Goldberg and A. V. Karzanov. Maximum skew-symmetric flows and matchings. Math. Program., 

Ser. A, 100:537-568, 2004. 

[54] A. V. Goldberg and R. Kennedy. Global price updates help. SIAM J. Discrete Mathematics, 10(4):551- 
572, 1997. 

[55] F. Hadlock. Finding a maximum cut of a planar graph in polynomial time. SIAM J. Comput, 
4(3):221-225, 1975. 

[56] M. Hagemann and O. Schenk. Weighted matchings for preconditioning symmetric indefinite linear 
systems. SIAM J. Sci. Comput, 28(2):403-420, 2006. 

[57] Yijie Han. Deterministic sorting in 0(n log log n) time and linear space. In Proceedings 34th ACM 
Symposium on Theory of Computing (STOC), pages 602-608. ACM Press, 2002. 

[58] S. Hanke. Zwei approximative Algorithmen fiir das Matchingproblem in gewichteten Graphen. PhD 
thesis, Humboldt-Universitat zu Berlin, 2004. 

[59] S. Hanke and S. Hougardy. New approximation algorithms for the weighted matching problem. Re- 
search Report No. 101010, Research Institute for Discrete Mathematics, University of Bonn, 2010. 

[60] N. Harvey. Algebraic algorithms for matching and matroid problems. SIAM J. Comput, 39(2):679-702, 
2009. 

[61] B. Hendrickson and R. Leland. The Chaco user's guide: Version 2.0. Technical Report SAND94-2344, 
Sandia National Laboratories, 1995. 



32 



[62] F. L. Hitchcock. The distribution of a product from several sources to numerous localities. J. Math. 

Physics, 20:224-230, 1941. 

[63] A. J. Hoffman and H. M. Markowitz. A note on shortest path, assignment, and transportation problems. 
Naval Research Logistics Quarterly, 10(l):375-379, 1963. 

[64] M. Holtgrewe, P. Sanders, and C. Schulz. Engineering a scalable high quality graph partitioner. In 
Proceedings 24th IEEE Symposium on Parallel and Distributed Processing (IPDPS), pages 1-12, 2010. 

[65] J. E. Hopcroft and R. M. Karp. An n^/^ algorithm for maximum matchings in bipartite graphs. SIAM 
J. Comput, 2:225-231, 1973. 

[66] C.-C. Huang and T. Kavitha. Efficient algorithms for maximum weight matchings in general graphs 
with small edge weights. In Proceedings ACM-SIAM 23rd Symposium on Discrete Algorithms (SODA), 
pages ??-??, 2012. 

[67] O. H. Ibarra and S. Moran. Deterministic and probabilistic algorithms for maximum bipartite matching 
via fast matrix multiplication. Info. Proc. Lett., 13(1):12-15, 1981. 

[68] P. Indyk. A near linear time constant factor approximation for Euclidean bichromatic matching (cost). 
In Proceedings 18th ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 39-42, 2007. 

[69] M. Iri, K. Murota, and S. Matsui. Heuristics for planar minimum-weight perfect matchings. Networks, 
13(l):67-92, 1983. 

[70] D. B. Johnson. Priority queues with update and finding minimum spanning trees. Info. Proc. Lett., 

4(3):53-57, 1975. 

[71] L. Kantorovitch. On the translocation of masses. Doklady Akad. Nauk SSSR, 37:199-201, 1942. 

[72] M.-Y. Kao, T.-K. Lam, W.-K. Sung, and H.-F. Ting. A decomposition theorem for maximum weight 
bipartite matchings. SIAM J. Comput., 31(l):18-26, 2001. 

[73] G. Karypis and V. Kumar. A fast and high quality multilevel scheme for partitioning irregular graphs. 
SIAM J. Sci. Comput, 20(l):359-392, 1998. 

[74] A. V. Karzanov. An exact estimate of an algorithm for finding a maximum flow, applied to the problem 

"on representatives" [in Russian]. Problems in Cybernetics, 5:66 70, 1973. Announced at the Seminar 
on Combinatorial Mathematics (Moscow, 1971). English translation available at the author's website. 

[75] A. V. Karzanov. Eflicient implementations of Edmonds' algorithms for flnding matchings with max- 
imum cardinality and maximum weight. In A. A. Pridman, editor. Studies in Discrete Optimization, 
pages 306-327. Nauka, Moscow, 1976. 

[76] H. W. Kuhn. The Hungarian method for the assignment problem. Naval Research Logistics Quarterly, 
2:83-97, 1955. 

[77] H. W. Kuhn. On combinatorial properties of matrices. George Washington University Logistics Pa- 
pers, 11:1-11, 1955. English translation of J. Egervary, Matrixok kombinatorius tulajdonsagairol, 
Matematikai es Fizikai Lapok 38, 16-28, 1931. 

[78] H. W. Kuhn. Variants of the Hungarian method for assignment problems. Naval Research Logistics 
Quarterly, 3:253-258 (1957), 1956. 

[79] H. W. Kuhn and W. J. Baumol. An approximate algorithm for the fixed-charges transportation 
problem. Naval Research Logistics Quarterly, 9:1-15, 1962. 



33 



[80] J. M. Kurtzberg. On approximation methods for the assignment problem. J. ACM, 9(4):419-439, 
1962. 

[81] M. K. Kwan. Graphic programming using odd or even points. Chinese Mathematics, 1:273-277, 1962. 

[82] E. Lawlcr. Combinatorial Optimization: Networks and Matroids. Holt, Rinehart & Winston, New 
York, 1976. 

[83] E. Leonard!, M. Mellia, F. Neri, and M. A. Marsan. Bounds on delays and queue lengths in input- 
queued cell switches. J. ACM, 50(4):520-550, 2003. 

[84] N. McKeown. The iSLIP scheduling algorithm for input-queued switches. IEEE/ACM Transactions 
on Networking, 7(2): 188-201, 1999. 

[85] N. McKeown, V. Anantharam, and J. C. Walrand. Achieving 100% throughput in an input-queued 
switch. In Proceedings 15th INFOCOM, pages 296-302, 1996. 

[86] S. Micali and V. V. Vazirani. An • algorithm for finding maximum matching in general 

graphs. In Proceedings 21st IEEE Symposium on Foundations of Computer Science (FOCS), pages 
17-27, 1980. 

[87] T. S. Motzkin. The assignment problem. Proceedings of Symposia in Applied Mathematics, VI, Nu- 
merical Analysis:lQ^12f), 1956. 

[88] M. Mucha and P. Sankowski. Maximum matchings via Gaussian elimination. In Proceedings 45th 
IEEE Symposium on Foundations of Computer Science (FOCS), pages 248-255, 2004. 

[89] J. Munkres. Algorithms for the assignment and transportation problems. J. Soc. Indust. Appl. Math., 
5:32-38, 1957. 

[90] F. OUivier. Looking for the order of a system of arbitrary ordinary differential equations. Appl. Algebra 
Eng. Commun. Comput., 20(l):7-32, 2009. English translation of: C. G. J. Jacobi, De investigando 
ordine systematis aequationum different ialum vulgarium cujuscunque, Borchardt Journal fiir die reine 
und angewandte Mathematik 65, 4, pp. 297-320, 1865, also reproduced in: C. G. J. Jacobi, Gesammelte 
Werke, Vol. 5, K. Weierstrass, Ed., Berlin, Bruck und Verlag von Georg Reimer, pp. 193-216, 1890. 

[91] M. Olschowka and A. Neumaier. A new pivoting strategy for Gaussian elimination. Linear Algebra 
and its Applications, 240:131-151, 1996. 

[92] J. B. Orlin and R. K. Ahuja. New scaling algorithms for the assignment and minimum mean cycle 
problems. Math. Program., 54:41-56, 1992. 

[93] G. I. Orlova and Y. G. Dorfman. Finding the maximum cut in a planar graph. Engineering Cybernetics, 
10:502 506, 1972. 

[94] F. Pellegrini. SCOTCH 5.1 user's guide. Technical report, LaBRI, 2008. 

[95] S. Pettie. Sensitivity analysis of minimum spanning trees in sub-inverse- Ackermann time. In Proceed- 
ings 16th Int'l Symposium on Algorithms and, Computation (ISAAC), pages 964 973, 2005. 

[96] S. Pettie and P. Sanders. A simpler linear time 2/3 — e approximation to maximum weight matching. 
Info. Proc. Lett., 91(6):271-276, 2004. 

[97] R. Preis. Linear time 1/2-approximation algorithm for maximum weighted matching in general graphs. 
In Proceedings 16th Symposium on Theoretical Aspects of Computer Science (STACS), LNCS 1563, 
pages 259-269, 1999. 



34 



[98] R. Preis and R. Dickmann. PARTY A software library for graph partitioning. In B. H. V. Topping, 
editor, Advances in Computational Mechanics with Parallel and Distributed Processing, pages 63-71. 
Civil-Comp Press, 1997. 

[99] E. M. Reingold and K. J. Supowit. Probabilistic analysis of divide-and-conquer heuristics for minimum 
weighted Euclidean matching. Networks, 13:49-66, 1983. 

[100] E. M. Reingold and R. E. Tarjan. On a greedy heuristic for complete matching. SIAM J. Comput., 
10(4):676-681, 1981. 

[101] J. Robinson. On the Hamiltonian game (a traveling salesman problem). Research Memorandum 303, 
RAND Corporation, 1949. 

[102] P. Sanders and C. Schulz. Distributed evolutionary graph partitioning. In Proceedings 14th Workshop 
on Algorithm Engineering and Experiments (ALENEX), pages 43-54, 2012. 

[103] P. Sankowski. Maximum weight bipartite matching in matrix multiplication time. Theoretical Com- 
puter Science, 410(44) :4480-4488, 2009. 

[104] O. Schenk, A. Wachter, and M. Hagemann. Matching-based preprocessing algorithms to the solution 
of saddle-point problems in large-scale nonconvex interior-point optimization. Comput. Optim. Appl., 

36(2-3):321 341, 2007. 

[105] D. Shah, P. Giaccone, and B. Prabhakar. Efficient randomized algorithms for input-queued switch 
scheduling. IEEE Micro, 22(1):10-18, 2002. 

[106] D. Shah and M. Kopikare. Delay bounds for the approximate maximum weight matching algorithm 
for input queued switches. In Proceedings 21st INFOCOM, pages 1024-1031, 2002. 

[107] R. L. Thorndike. The problem of classification of personnel. Psychometrika, 15:215-235, 1950. 

[108] M. Thorup. Undirected single-source shortest paths with positive integer weights in linear time. 
J. ACM, 46(3):362-394, 1999. 

[109] M. Thorup. Integer priority queues with decrease key in constant time and the single source shortest 
paths problem. In Proceedings 35th ACM Symposium on Theory of Computing (STOC), pages 149-158, 
2003. 

[110] M. Thorup. Equivalence between priority queues and sorting. J. ACM, 54(6), 2007. 

[Ill] N. Tomizawa. On some techniques useful for solution of transportation network problems. Networks, 
1(2):173-194, 1971. 

[112] K. R. Varadarajan and P. K. Agarwal. Approximation algorithms for bipartite and non-bipartite 
matching in the plane. In Proceedings 10th ACM-SIAM Symposium on Discrete Algorithms (SODA), 
pages 805-814, 1999. 



[113] V. V. Vazirani. A theory of alternating paths and blossoms for proving correctness of the 0{\/VE) 
general graph maximum matching algorithm. ComMnatorica, 14(1) :71 109, 1994. 

[114] D. E. D. Vinkemeier and S. Hougardy. A linear-time approximation algorithm for weighted matchings 
in graphs. ACM Trans, on Algorithms, 1(1):107-122, 2005. 

[115] J. von Neumann. A certain zero-sum two-person game equivalent to the optimal assignment problem. 
In H. W. Kuhn and A. W. Tucker, editors. Contributions to the Theory of Games, volume II, pages 
5-12. Princeton University Press, 1953. 



35 



[116] C. Walshaw and M. Cross. JOSTLE: Parallel multilevel graph-partitioning software an overview. In 

Mesh Partitioning Techniques and Domain Decomposition Techniques, pages 27-58, 2007. 

[117] C. Witzgall and C. T. Zahn, Jr. Modification of Edmonds' maximum matching algorithm. J. Res. 
Nat. Bur. Standards Sect. B, 69B:91-98, 1965. 



36 



